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EXECUTIVE  SUMMARY 


This  document  is  a  set  of  notes  intended  to  help  a  new  Ph.D.  in  physics  or 
chemistry  (or  person  of  equivalent  training)  become  familiar  with  the  specialized  literature 
of  the  applied  fields  of  atmospheric  nuclear  weapons  effects,  missile  reentry,  upper 
atmospheric  physics,  or  high-energy  laser  propagation  through  the  atmosphere.  The 
document  is  not  a  substitute  for  a  course  in  high-temperature  gas  dynamics.  It  presents  a 
simple  overview  of  some  key  physical  factors,  and  gives  references  both  to  relevant 
textbooks  and  to  the  "gray"  literature  specific  to  the  field.. 

Part  I  consists  of  a  general  introduction  (Chapter  1),  a  discussion  of  equilibrium  air 
at  temperatures  to  approximately  10,000  K  (Chapter  2),  a  discussion  of  the  radiation  from 
high-temperature  air  (Chapter  3),  and— as  a  bridge  to  Part  H,  to  be  published  later,  which 
will  address  some  applications— Chapter  4,  which  gives  a  brief  discussion  of  some  specific 
heating  mechanisms. 

Chapter  2  points  out  that  air  starts  to  dissociate  at  a  temperature  of  2000  K  to 
4000  K  (depending  on  the  pressure),  even  though  the  mean  thermal  energy  of  the 
molecules  at  these  temperatures  is  much  less  than  the  dissociation  energy.  The  emphasis  in 
the  discussion  is  on  equilibrium  processes,  with  an  indication  of  where  detailed 
considerations  of  non-equilibrium  such  as  chemical  kinetics  become  important.  While 
many  current  applications  deal  with  non-equilibrium  situations,  yet  those  generally  start 
from  equilibrium  which  provides  a  simple  and  rigorous  limit.  It  is  pointed  out  that  at  the 
higher  temperatures  T,  radiative  energy  transfer  (which  goes  as  T^  or  faster)  tends  to 
predominate  over  convective  heat  transfer  (which  goes  as  grad  T).  Chapter  3  treats 
radiative  problems,  addressing  both  radiative  transfer  and  quantitative  molecular 
spectroscopy,  and  finishes  by  comparing  radiative  heat  transfer  by  dissociating  air  and  by 
other  molecules  with  that  from  a  blackbody  in  the  same  temperature  range.  Chapter  4 
addresses  the  thermal  excitation  of  air  by  shock  and  other  aerodynamic  heating  (in 
spacecraft  reentry),  by  X-ray  heating  (in  the  nuclear  fireball),  and  by  the  absorption  of  solar 
ultraviolet  radiation  (in  the  earth's  upper  atmosphere). 
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Part  n,  Applications,  will  be  prepared  later.  It  will  expand  on  Chapter  4  with  more 
extensive  discussions  of  reentry  physics,  atmospheric  nuclear  weapons  effects,  and 
electromagnetic  wave  propagation  within  the  atmosphere. 


S-2 


1.  INTRODUCTION 


A  variety  of  technical  applications  require  an  understanding  of  the  physics  of  high- 
temperature  air,  where  "high  temperature"  is  taken  to  mean  1000  K  to  10,000  K  .  Some 
applications  are  to  missile  or  space  vehicle  entry  into  the  earth's  atmosphere,  and  to  nuclear 
explosions  within  the  atmosphere  which  lead  to  the  formation  of  a  "fireball."  Related 
applications  deal  with  rocket  propulsion  and  with  the  propagation  of  high-energy  laser 
beams  in  the  atmosphere. 

When  air  is  heated  to  2000  K  to  4000  K  (depending  on  pressure)  it  begins  to 
dissociate;  by  10,000  K  it  is  not  only  fully  dissociated  but  at  least  singly  ionized.  To  break 
all  these  chemical  bonds  requires  a  great  deal  of  energy.  Thus  the  physical  and  chemical 
properties  of  high-temperature  air  are  quite  different  from  those  of  room  temperature  air 
because  of  the  presence  of  reactive  molecules  (such  as  oxides  of  nitrogen),  atoms,  and 
(atomic  and  molecular)  ions. 

The  general  training  of  a  Ph.D.  physicist  or  chemist  typically  does  not  Include  the 
characteristics  of  air  in  the  1000  K  to  10,000  K  temperature  range.  These  properties  were 
elucidated  largely  between  1950  and  1965,  frequently  in  an  applied  and  often  classified 
context  in  the  "gray  literature"  of  technical  reports.  Someone  who  begins  working  in  the 
field  of  reentry  physics  now,  in  1989,  may  be  referred  to  computer  codes  which  were 
generated  originally  as  much  as  25  years  ago,  have  been  modified  extensively  and  are 
probably  documented  and  referenced  inadequately.  There  have  been  very  few  atmospheric 
nuclear  tests  since  1962,  and  the  literature  consists  of  successively  more  sophisticated 
numerical  analyses  of  these  very  limited  test  data.  In  either  case,  the  physical  basis  of 
current  numerical  models  is  often  not  readily  accessible  to  the  new  user. 

In  much  of  the  current  literature,  the  chemical  composition  of  the  relevant  gases  is 
based  on  the  solution  of  systems  of  coupled  chemical  reactions  rather  than  on  the 
assumption  of  Local  Thermodynamic  Equilibrium  (LTE).  Such  a  non-equilibrium 
description  can  be  more  accurate  than  LTE  for  times  so  short  that  equilibrium  is  not 
attained,  provided  the  reaction  set  used  is  sufficiently  large,  appropriate,  and  correct,  which 
is  a  non-trivial  condition.  However,  as  an  introduction  to  the  field  it  is  convenient  to 
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discuss  the  physics  largely  in  terms  of  an  LTE  model  which  is  intuitively  simple  and 
provides  a  starting  point  for  the  more  detailed  non-LTE  models. 

This  document  is  the  first  part  of  a  two-part  series.  Discussed  here  are  the 
equilibrium  composition  and  radiative  properties  of  high-temperattue  air  (respectively,  in 
Chapters  2  and  3),  indicating  where  non-equilibrium  (kinetic)  considerations  become 
important.  Chapter  4  is  a  bridge  to  a  planned  Part  II  of  this  document  which  will  be 
published  later.  It  presents  some  examples  of  how  the  heated  air  discussed  here  can 
actually  be  produced:  by  aerodynamic  heating  in  reentry,  where  the  kinetic  energy  of 
hypersonic  motion  is  transformed  into  random  molecular  motion  or  heat,  by  the  absorption 
of  X-rays  in  the  atmsaophere  that  produces  nuclear  fireballs,  or  by  the  absorption  of  solar 
ultraviolet  radiation  in  the  earth's  upper  atmosphere,  which  produces  the  ionosphere/ 
thermosphere. 

Part  n  will  deal  in  more  detail  with  missile  reentry  (Chapter  5),  with  nuclear 
phenomenology  (Chapter  6),  and  with  selected  topics  in  the  propagation  of  visible  and 
infrared  radiation  from  high-temperature  sources  through  the  atmosphere  (Chapter  7). 

This  document  is  a  set  of  notes  intended  to  help  a  new  Ph.D.  in  physics  or 
chemistry  (or  person  of  equivalent  training)  become  familiar  with  the  specialized  liteiature 
of  an  applied  field  such  as  atmospheric  nuclear  weapons  effects,  missile  reentry,  upper 
atmospheric  physics,  or  high-energy  laser  propagation  through  the  atmosphere.  It  presents 
a  simple  overview  of  some  key  physical  factors,  and  gives  references  both  to  relevant 
textbooks  and  to  the  "gray"  literature  specific  to  the  field. 

While  this  document  has  been  read  critically  by  a  number  of  people— see  the 
acknowledgments— obscurities  and  inaccuracies  probably  remain.  The  reader  should  be 
aware  of  these  limitations  in  using  the  material,  and  is  invited  to  give  comments  and 
criticisms  to  the  author  to  help  make  the  material  more  useful. 
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2.  COMPOSITION  OF  HIGH-TEMPERATURE  AIR 


2.1  INTRODUCTION 

2.1.1  Basic  Data  for  Air 

When  air  at  normal  pressures  (say,  10“^  to  1  atmospheres)  is  heated  from  room 
temperature  to  perhaps  10,000  K,  most  of  the  energy  absorbed  goes  into  dissociating 
oxygen  molecules  rather  than  into  increasing  the  kinetic  energy  of  the  O2  and  N2 
molecules.^  This  fact  is  a  key  motivation  in  preparing  the  present  set  of  notes,  because  the 
physics  in  the  regime  discussed  involves  the  excitation  of  various  internal  atomic  and 
molecular  degrees  of  freedom,  in  particular  molecular  dissociation. 

Note  that  for  T  =  10,000  K,  kT/e  =  0.86  eV  (where  k  =  Boltzmann's  constant  and 
e  =  charge  of  the  electron).  The  dissociation  energy  of  the  O2  molecule  is  D  =  5.1 1  eV 
which  corresponds  to  a  characteristic  temperature  0d  =  eP/k  =  59300  K,  but  at  normal 
pressures  O2  is  dissociated  already  at  temperatures  of  4000-5000  K,  in  essence  because  the 
phase  space  volume  of  O  +  O  is  much  greater  than  that  of  02^ 

Table  2-1  lists  the  composition  of  dry  air  at  sea  level.  Notice  that  the  chemically 
and  radiatively  reactive  and  therefore  most  interesting  species— CO2,  CO,  O3,  NO,  and 
NO2  (also  CH4,  NH3,  and  SO2  which  are  not  discussed  here)— are  present  in  small  and 
typically  quite  variable  concentrations. 

Table  2-2  gives  some  properties  of  air  up  to  120  km  altitude,  where  most  molecules 
other  than  N2  are  significantly  dissociated.^  This  table  lists  the  temperature,  density,  gas- 
kinetic  mean  free  path  (MFP)  as  well  as  the  mean  mixing  ratios  for  the  energetically, 
radiatively,  and  chemically  important  minor  species  H2O  and  O3.  Note  the  following: 


'  See  Section  4.2,  in  particular  Table  4-1. 

^  See  the  discussion  of  Section  2.2.2. 

3  Note  that  this  is  due  to  solar  photodissociation,  not  to  thermal  dissociation.  Section  4.4  gives  a  brief 
discussion. 


Table  2-1.  Composition  of  Dry  Air  at  1  Atmosphere  Pressure 
(From  U.S.  Standard  Atmosphere  1976,  pp.  3,  33,  38,  44) 


Molecule 

Molecular  Weight 

Volume  Fraction 

N2 

28 

0.781 

02 

32 

0.209 

Ar 

39.9 

0.0093 

CO2 

44 

0.00032  (increasing  somewhat  after  1980) 

Ne 

20.2 

0.000018 

He 

4.0 

0.0000052 

Kr 

83.8 

0.0000011 

CH4 

16.0 

1500  X  10-® 

H2 

2.02 

500  X 10-9 

CO 

28 

190x10-9 

O3 

48 

40  X  10“9  (highly  variable~see  p.  39,  op.  cit.) 

NH3 

17 

4x10-9 

SO2 

64 

1  X  10-9 

NO2 

46 

1  X  10-9 

NO 

30 

1 

0.5  X  10-9 
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Table  2-2.  The  Atmosphere  to  High  Altitudes 


_  I 


■BB 

■oSBi 

mSSSm 

Comments 

120 

360 

2.22E-8 

3.23 

0.20 

0.0005 

O2  dissoc.  (A) 

110 

240 

9.71  E-8 

0.815 

0.28 

0.05 

O2  dissoc.  (A) 

100 

195 

5.6E-7 

0.163 

0.40 

0.4 

90 

187 

3.42E-6 

0.0256 

0.85 

0.7 

80 

199 

1.85E-5 

4.07E-3 

2.1 

^B 

T  min.  (B) 

70 

220 

8.28E-5 

9.28E-4 

3.5 

■■ 

60 

247 

3.1  E-4 

5.52E-4 

4.8 

1.1 

50 

271 

1.03E-3 

7.91  E-5 

5.2 

3.1 

T  max.  (C) 

40 

250 

4.00E-3 

2.03E-5 

5.0 

7.3 

30 

227 

0.0184 

4.41  E-6 

4.7 

6.6 

20 

217 

0.0889 

9.14E-7 

3.9 

2.6 

Low  H2O  (D) 

10 

223 

0.414 

1 .96E-7 

70(E) 

0.13 

T  min.  (F) 

0 

288 

1.225  I 

! 

6.63E-8 

7800  (E) 

0.027 

Notes: 

*  ppmV  >  parts  per  million  by  volume 

*  Temperature,  density,  mean  free  path  (MFP)  from  USSA-62--see  USAF,  1 965,  p.  2-1 9. 


*  H2O,  O3  from  AFGL  atmospheric  profiles-see  Anderson  et  al.,  1 986. 

(A)  O2  begins  to  dissociate  above  ~  80  km.  By  110-120  km,  the  number  density  of  atomic  O  equals 
that  of  02-  Ozone,  water  vapor,  and  carbon  dioxide  all  dissociate  at  somewhat  lower  altitudes. 

(B)  Temperature  minimum  (mesopause)  near  85  km. 

(C)  Temperature  maximum  (stratopause)  near  50  km. 

(D)  The  atmosphere  is  very  dry  above  the  trr^pause  {-  1 1  km). 

(E)  The  troposphere  (below  -•  1 1  km)  is  quite  moist,  with  a  highly  variable  humidity. 

(F)  Temperature  minimum  (tropopause)  near  1 1  km  at  mid-latitudes. 
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•  The  mean  free  path  [defined  in  Section  2.5— see  Eq.  (2.13)]  is  7  x  10^  cm  = 
700  A  at  sea  level,  but  16  cm  at  100  km  altitude.  The  MFP  is  proportional  to 
the  inverse  of  the  number  density,  and  normal  (~  1  m  )  reentry  vehicles  will 
begin  to  interact  with  the  atmosphere  at  altitudes  below  1 10  to  100  km. 

•  The  lower  atmosphere  (troposphere,  below  the  temperature  minimum  of 
tropopause  near  1 1  km  at  rmd-latitudes)  is  quite  moist,  with  highly  variable 
humidity.  The  minimum  temperature  of  ~  217  K  at  the  tropopause  provides  a 
"cold  trap"  which  reduces  the  mixing  ratio  above  the  tropopause.'^ 

•  The  ozone  maximum  occurs  in  the  stratosphere  as  O2  is  photodissociated  at 
higher  altitudes  while  one  of  the  0-atoms  produced  in  this  way  bonds  to 
another  O2  molecule  to  produce  ozone,  which  is  then  photodissociated  and 
reformed  repeatedly. 

Table  2-3  lists  dissociation  and  ionization  potentials  and  some  other  information  for 
important  air  molecules  and  atoms.  This  basic  information  is  of  importance  in  the  present 
context.  Note  in  particular 

•  The  N2  molecule  has  a  significantly  higher  dissociation  eneregy  (9.76  eV)  than 
any  of  the  other  molecules,  which  explains  the  persistence  of  this  chemical 
species. 

•  The  O2  molecule  has  the  lowest  ionization  potential  of  any  of  the  significant 
species,  which  explains  why  this  is  often  the  predominant  positive  ion.^ 

•  Vibration  frequencies  are  normally  in  the  range  l(X)0-20()0  cm~^  which  is 
why  the  molecular  vibration  spectra  occur  in  the  IR,  as  the  corresponding 
wavelengths  are  10-  5^im. 

•  Molecular  rotation  constants  Be  are  normally  0.5-2  cm~l  except  for  hydrogen- 
containing  species  (H2O)  which  have  a  much  lower  rotational  constant  (and 
larger  vibrational  frequencies  cofe). 

Some  information  about  the  quantum  mechanical  strucmre  of  atoms  and  air 
molecules  is  listed  in  Appendix  A.  For  more  detail,  refer  to  the  bibliography. 


Some  numbers  are  in  order.  The  saturated  water  vapor  pressure  over  ice  at  217  K  is  0.018  mb  which 
corresponds  to  a  maximum  mixing  ratio  at  1 1  km  of  80  ppm  V.  By  contrast,  the  saturated  water  vapor 
pressure  over  water  at  288  K  is  17  mb  which  gives  a  maximum  sea  level  mixing  ratio  of 
17,000  ppmV  (cf.  Smithsonian  Met.  Tables,  1949.  p.  234  0- 

^  NO  and  NO2  have  lower  ionization  potentials,  but  these  molecules  are  typically  present  in  relatively 
low  concentration. 
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Table  2-3.  Some  Properties  of  Principal  Air  Atoms  and  Molecules 


Species 

D 

1 

(Oe 

Be 

02(x3lg-) 

5.116 

12.07 

1580 

1.446 

N2{X1V) 

9.76 

15.58 

2359 

1.998 

NO  (X2nr) 

6.497 

9.26 

1904 

1.672 

C02(X1V) 

5.45  (O-CO) 

13.8 

1388,667,2349 

0.39 

03(X1Ai) 

1.04  (O-O2) 

12.5 

1103,701,1042 

3.55,0.445,0.39 

N02(X2Ai) 

3.11  (O-NO) 

9.8 

1320,750,1618 

8.0,0.43,0.41 

4.50  (N-O2) 

H20(X1Ai) 

5.12  (H-OH) 

12.6 

3657,1595,3756 

27.9,14.5,9.29 

0(3P) 

13.62 

N(4S) 

14.53 

Ar(1S) 

15.76 

Notes: 


0  (eV)  -  dissociation  energy. 

I  (eV)  -  ionization  potential 
(Oe  (cm“^ )  »  vibrational  frequency. 
Be  (cm“^ )  »  rotational  constant. 


For  notation  on  electronic  states  see  Appendix  A. 


2.1.2  On  Thermodynamic  Equilibrium 

Air  is  made  up  of  a  number  of  chemically  reacting  species,  and  thus  its  composition 
under  given  conditions  is  determined  by  the  detailed  balance  of  different  chemical  reaction 
processes.  For  a  complex  system,  a  large  number  of  reactions  may  have  to  be  considered. 
While  the  necessary  computations  can,  in  general,  be  carried  out,  given  present-day 
computer  power,  one  is  sometimes  uncertain  whether  all  important  reactions  have  been 
included,  and  whether  their  rate  constants  are  represented  correctly. 

If  we  ask  for  the  state  of  a  gas  after  a  sufficiently  long  time,  so  that  the  detailed 
collisional  and  radiative  energy  transfer  processes  have  had  an  opportunity  to  establish  a 
balance,  then  many  properties  of  the  gas  are  characterized  simply  as  functions  of 
temperature  and  pressure,  or  of  temperature  and  density,  independent  of  details  of  reaction 
mechanisms  and  rate  constants.  This  is  the  very  useful  concept  of  Thermodynamic 
Equilibrium,  which  is  a  very  useful  limit  for  high-temperature  air  in  the  applications 
considered  here  (reentry  and  nuclear  fireballs),  and  it  is  a  good  starting  point  for  more 
complicated  applications. 

Even  if  the  Umit  of  thermodynamic  equilibrium  does  not  hold  universally,  it  may 
still  apply  for  some  specific  processes,  such  as  the  excitation  of  rotational  levels  of  a  given 
molecule  or  the  dissociation  of  one  or  another  molecule,  and  thus  we  can  use  the  concept  of 
Quasi- Equilibrium  for  the  particular  species  to  which  it  applies. 

However,  if  we  ask  for  the  changing  properties  of  a  medium  at  times  so  early  that 
thermodynamic  equilibrium  has  not  yet  been  achieved,  it  is  necessary  to  consider  the  details 
of  the  Kinetic  Regime  (involving  collisional  and  radiative  processes)  that  underlies  the  state 
of  the  particular  medium.  While  it  may  be  necessary  to  go  to  this  limit  in  an  actual 
application,  this  will  not  be  discussed  here,  since  for  our  purposes  it  is  better  to  begin  with 
the  simpler  case  of  (near)  equilibrium. 

2.1.3  Outline  of  the  Chapter 

In  Section  2.2  we  discuss  the  equilibrium  dissociation  of  O2  and  the  equilibrium 
ionization  of  atomic  oxygen,  whose  ionization  potential  is  13.62  eV.  It  is  largely  ionized 
by  12,000  K  to  15,000  K  where  kT  =  1  -  1.3  eV.  This  analysis  demonstrates  how  the 
degrees  of  dissociation  (and  ionization)  vary  with  the  dissociation  (and  ionization)  energy, 
temperature,  and  pressure  or  density.  Because  only  a  single  chemical  process  is 
considered,  the  analysis  is  simple.  In  Section  2.3  some  results  are  presented  on  the 
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equilibrium  composition  of  air  at  elevated  temperatures  obtained  by  minimizing  the 
(Helmholtz)  free  energy  numerically.  The  basic  results  which  are  still  in  general  use  are 
more  than  20  years  old— newer  values  of  thermodynamic  input  parameters  provide  slight 
improvements^  which  may  be  important  for  particular  applications.  Because  the 
dissociation  energy  is  so  important  in  the  regime  being  considered  here.  Section  2.4 
discusses  some  of  its  effects  on  thermodynamic  properties. 

The  discussion  of  Sections  2.2  through  2.4  considers  thermodynamic  equilibrium, 
but  this  is  not  uniformly  applicable  for  chemically  and  radiatively  important  species.  The 
concentration  of  water  vapor—while  small— is  highly  variable,  and  the  important  species  NO 
is  metastable:  it  is  formed  in  amounts  of  several  percent  at  temperatures  above  2000  K,  and 
at  lower  temperatures  its  chemical  disappearance  rate  is  so  slow  that  it  tends  to  remain 
without  decomposing  ("freezes-in")  as  the  air  cools.^  Thus,  Section  2.5  treats  some  non¬ 
equilibrium  processes,  particularly: 

•  When  kinetics  are  important  (low  pressure,  uncommon  species,  short  time). 

•  Some  peculiarities  of  the  NO  molecule,  which  is  an  important  radiating 
species,  has  a  low  ionization  potential,  and  is  chemically  more  reactive  than  N2 
or  O2. 

•  Some  qualitative  factors  on  radiating  species. 

Radiative  heat  transfer  tends  to  dominate  over  convective  heat  transfer  for  gases  at 
very  high  temperatures,  and  thus  Section  2.6,  which  sketches  some  heat  transfer  issues  of 
importance  in  high-temperature  gases,  serves  as  an  introduction  to  the  analysis  of  radiative 
processes  in  Chapter  3. 


^  The  basic  work  on  the  thermodynamic  properties  of  air  species  was  done  in  the  1950s  and  1960s;  for 
more  recent  work,  see,  e.g..  Moss  and  Scott,  1985. 

^  The  term  "freeze-in"  is  used  in  the  specialized  literature  of  this  field  to  indicate  that,  as  the  temperature 
of  a  gas  parcel  falls,  its  chemical  composition  no  longer  follows  thermodynamic  equilibrium,  but 
rather  some  high-temperature  species  (such  as  NO,  O,  or  OH)  may  survive  for  a  relatively  long  time  at 
lower  temperatures  because  the  decomposition  processes  for  these  high-temperature  species  may 
become  very  slow.  See,  e.g.,  Zeldovich  and  Raizer,  1967,  p.  564f. 
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2.2  EQUILIBRIUM  DISSOCIATION  AND  IONIZATION  FOR  A 

SINGLE  SPECIES 

2.2.1  Introduction 

Here  we  present  the  simple  expressions  and  numerical  results  for  the  dissociation  of 
the  oxygen  molecule  and  for  the  ionization  of  the  oxygen  atom.  The  results  are  idealized  in 
that  the  excited  states  of  O  and  O2  listed  in  Appendix  A  are  not  considered.  These 
expressions  indicate  in  parametric  form  how  the  dissociation  and  ionization  depend  on 
temperature  and  density  (or  pressure),  as  well  as  on  the  effective  dissociation  or  ionization 
energy.  They  shows  that  at  reasonable  pressures  ( less  than  1  atmosphere)  the  O2  molecule 
of  dissociation  energy  D  =  5. 1  eV  is  already  significantly  dissociated  at  a  temperature  of 
5000  K  which  is  less  than  10  percent  of  the  effective  dissociation  temperature  ©d  =  D/k  = 
59,300  K;  also  the  0-atom  of  ionization  potential  I  =  13.6  eV  is  significantly  ionized  at  a 
temperature  of  12,000  K  which  is  less  than  10  percent  of  the  effective  ionization 
temperature  ©i  =  I/k  =  158,000  K. 

2.2.2  Dissociation  of  Molecular  Oxygen 

Consider  the  diatomic  molecule  O2  with  dissociation  energy  D  which  is  made  up  of 
two  O  atoms.  At  temperature  T,  if  a  total  number  Ntot  of  molecules  occupying  volume  V 
are  made  up  of  N02  molecules  and  Nq  O-atoms,  then  in  local  thermodynamic  equilibrium 
the  respective  numbers  No  and  N02  satisfy  the  following  two  equations: 

N02+N0/2  =  Ntot  =  constant  (2.1) 

which  is  the  equation  of  continuity ,  and  the  chemical  equilibrium  condition 

No2/No2  =  KdissCT)  =  [fo^T)  /fo2(T)]  exp(-D/kT)  (2.2) 

where  K<iiss(T)  is  the  equilibrium  constant  for  the  dissociation  reaction 

O2  <=>  20  (2.3) 

and  the  fs  are  partition  functions  for  the  respective  atoms  and  molecules  which  are 
proportional  to  volume  V  The  partition  function  fj(T)  for  a  system  J  having  states  i  = 
1,2,...  of  energy  Ej  is  defined  as  the  sum  over  all  states  i  of  the  factor  exp(-Ei/kT) 

fj(T)  =  Lj  exp(-Ei/kT)  .  (2.4) 
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<3 —  phi(1) 

♦  phi{1e-3) 

■■ -  phi(1e-6) 

phi(1e-9) 


Temp.(K) 


Figure  2>1.  Thermal  Equilibrium  Dissociation  of  Molecular  Oxygen  O2. 
as  a  Function  of  Temperature  T  and  of  Number  Density  Ratio  n/no. 

The  fractional  dissociation  phi  ((t>diss  of  Eq.(2.5))  Is  expressed  as  a  function  of 
the  number  density  ratio,  n/no,  where  no  =  2.687x10''^  cm*^. 
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The  partition  function  can  be  written  to  a  good  approximation  as  the  product  of 
translational,  rotational,  vibrational  and  electronic  terms;  in  Appendix  B  we  give  the 
partition  function  for  the  dissociation  of  a  diatomic  molecule  and  the  numerical  parameters 
for  the  dissociation  of  the  oxygen  molecule. 

In  Fig.  2-1  we  show  the  fractional  dissociation  of  oxygen,  defined  as: 

<l>diss  =  No/[No  +  No2]  (2.5) 

for  various  values  of  the  number  density  ratio  n/no,  where  no  =  2.687  x  10^^  particles/cm^, 
which  corresponds  to  a  pressure  of  1  atmosphere  at  a  temperature  of  273  K  for  an  ideal 
gas.® 

The  dissociation  energy  D  =  5.1 1  eV  can  be  expressed  as  an  effective  dissociation 
temperature  ©d  =  D/k  =  59,300  K.  Note  that  the  temperature  dependence  of  the  fractional 
dissociation  <|)diss  is  much  stronger  than  the  pressure  (or  specific  volume  or  density) 
dependence,  and  that  at  n/no  >  10~3*  which  correspond  to  a  pressure  of  less  than  0.015 
atmospheres,  the  O2  molecule  is  already  significantly  dissociated  at  4000  K.^ 

2.2.3  Ionization  of  Atomic  Oxygen 

For  definiteness,  consider  the  single  ionization  of  the  oxygen  atom: 

O  +  Energy  (I  =  13.6  eV)  =>  0+  +  e"  .  (2.6) 

We  begin  with  Najt  atoms  in  a  volume  V;  if  in  equilibrium  at  temperature  T  there  are  Na 
atoms,  Ni  ions  and  Ne  electrons,  the  condition  of  charge  neutrality  requires  that 

Ni  =  Ne  (2.7) 

and 

Na  +  Ni  =  Ntot  •  (2.8) 


®  It  is  appropriate  to  take  V  =  22.4  liters/mole  and  Ntot  =  (n/n©)  NAvogadro-  This  would  give  a  pressure 
(T/273)  (n/no)  (1  -  ')>diss)  at  a  temperature  T. 

^  The  reason  for  this  strong  temperature  as  against  density  dependence  can  be  seen  from  the  form  of  the 
equation  derived  for  Nq  from  Equations  (2.1)  and  (2.2),  which  is  No^/(Ntot-No)  =  G  V  expf-eo/T), 
where  G  has  cmly  a  weak  temperature  dependence  and  no  {xessuie  dependence). 
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The  equilibrium  constant  Ki(T)  for  the  ionization  reaction  (2.6)  is  given  by  the  relation 

Ni.Ne/NA  =  Ki(T)  =  [fi(T)  fe(T)/fA  (T)]  exp  (-I/kT)  ,  (2.9) 

where  the  f  s  are  partition  functions.  For  our  idealized  system,  Eq.  (2.9)  becomes 

Nc^/Na  =  V(27ckm‘/h2)3/2T3/2[2g(r)/g(A)]exp(-I/kT)  .  (2.10) 

Here  m'  is  the  reduced  mass  of  the  electron  in  presence  of  the  ion;  it  is  essentially  equal  to 
the  electronic  mass.  The  factor  2  associated  with  the  statistical  weights  g  comes  from  the 
two  spin  states  of  the  free  electron,  and  the  statistical  weights  g(I)  and  g(A)  are  respectively 
g[0‘'’(^S)]  =  4  and  g[0(3p)]  =  3.  The  rest  of  the  partition  function  comes  from  the  phase 
space  volume  term  associated  with  the  translational  partition  function  analogous  to 
Eq.  (B.la)  of  Appendix  B. 

We  define  the  degree  of  ionization  as 

4>ion  =  Ne/(NA  +  Ni)  ,  (2.11) 

where  Na  and  Ng  are  determined  by  solving  Eqs.  (2.8)  and  (2.10).  The  resulting  values 
for  <j)ion  as  a  function  of  temperature  and  density  are  plotted  in  Fig.  2-2  where  again  Ntot  = 
(n/no)  NAvogadro-  Note  that  once  again  the  temperature  dependence  of  (l)ion  is  much 
stronger  than  the  density  dependence,  and  that  at  temperatures  above  ©i/10=  16,000  K  the 
atomic  gas  is  largely  ionized  at  pressures  below  1  atmosphere. 

2.3  EQUILIBRIUM  COMPOSITION  OF  (DRY)  AIR  AS  FUNCTION  OF 
TEMPERATURE 

Air  is  made  up  of  O2  of  dissociation  energy  5.12  eV,  which  is  largely  dissociated  at 
"normal"  pressures  by  4,000  K  to  6,000  K,  and  of  N2,  whose  dissociation  energy  of 
D  =  9.76  eV  corresponds  to  ©d  =  cD/k  =  113,400  K,  and  which  does  not  begin  to 
dissociate  at  temperatures  below  ~  8000  K.  Also,  some  NO  is  formed  at  temperatures 
above  2000  K,  and  the  small  fraction  of  CO2  and  H2O  will  dissociate  in  the  range  1 ,000  K 
to  4,000  K.  To  calculate  the  composition  of  equilibrium  air  under  conditions  when  the 
chemical  composition  changes,  one  either  solves  a  set  of  equilibrium  conditions  analogous 
to  Eq.  (2.2)  or  (2.10)  for  each  of  the  reactions  considered,  or  equivalently  minimizes  the 
free  energy  at  constant  pressure  (using  the  Gibbs  free  energy)  or  at  constant  density  (using 
the  Helmholtz  free  energy).  This  is  done  with  a  computer  program,  using  as  input  the 
atomic  composition  and  the  thermodynamic  parameters  for  the  various  species,  some  of 
which  are  given  in  Table  2-3.  The  following  points  should  be  noted: 


2-11 


phi(l) 

phi(1e-3) 

phi(1e-6) 

phi(1e-9) 


Figure  2-2.  Thermal  Equilibrium  Ionization  of  Atomic  Oxygen  as  a  Function  of 
Temperature  T  and  of  Number  Density  Ratio  n/no-  The  fractional  Ionization  of 
phi  ((t)|on  of  Eq.(2.1l))  Is  expressed  as  a  function  of  the  number  density  ratio, 

n/no,  where  no  =  2.687x10^®  cm’®. 
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1 .  Most  calculations  have  been  done  for  dry  air  because  the  atmospheric  water 
vapor  concentration  is  small  and  highly  variable. 

2.  For  a  review  of  thermodynamic  calculations  in  the  1950s  and  early  1960s,  see 
Hochstim,  1964  (which  is  reprinted  as  an  appendix  to  Hilsenrath  and  Klein, 
1965). 

3.  The  "classic"  calculations  of  the  early  1960s  are  still  generally  accepted,  but 
see,  e.g..  Moss  and  Scott,  1985,  for  references  to  current  work  in  this  field. 

Some  representative  results  for  the  composition  of  high-temperature  dry  air  (from 
Hilsenrath  and  Klein,  1965)  are  shown  in  Fig.  2-3.  These  calculations  are  done  at  constant 
density,  i.e.,  with  p/po  =  10“",  where  po  =  1.293  kg/m^,  which  corresponds  to  STP.^® 
The  results  are  inevitably  complex,  but  some  general  conclusions  that  can  be  drawn  are  the 
following: 

•  The  air  triatomics  dissociate  below  2,000  K  to  3,000  K,  depending  on  density 
or  pressure. 

•  O2  (D  =  5.12  eV)  dissociates  in  the  3,000  K  to  6,000  K  temperature  range, 
while  N2  (D  =  9.76  eV)  dissociates  in  the  9,000  K  to  15,0(X)  K  range. 

•  There  is  little  ionization  below  6,000  K  where  NO+(which  has  the  lowest 
ionization  potential  of  any  of  the  air  species,  I  =  9.25  eV)  is  the  principal 
positive  ion;  above  9,000  K  there  are  very  few  molecules  left,  while  above 
12,000  K,  air  is  largely  ionized  [with  0+  (I  =  12.1  eV)  as  principal  positive 
ion]. 

•  Metastable  molecules  or  "free  radicals"  such  as  NO,  CO,  CN,  OH,  etc.,  which 
are  both  chemically  and  radiatively  active,  form  above  room  temperature  but 
normally  decompose  above  3,000  K  to  5,000  K. 

2.4  EFFECTS  OF  DISSOCIATION  ENERGY  ON  THERMODYNAMIC 

PROPERTIES 

The  characteristic  behavior  of  air  in  the  temperature  range  from  1,000  K  to 
10,000  K  which  is  under  discussion  here  is  principally  due  to  the  fact  that  a  large  amount 
of  energy  is  required  to  dissociate  the  O2  and  N2  molecules.  We  can  write  the  change  in 
internal  energy  AE  as 

AE  = -P  A V S  AT -f- p  An  ,  (2.8) 


^  ®  Po  corresponds  to  Hq  of  the  previous  Section. 
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Figure  2-3.  Equilibrium  Composition  of  Dry  Air  as  Function  of  Temperature 
and  Density.  Source:  Hilsenrath  and  Klein,  1965. 
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Figure  2-3.  (continued) 
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Figure  2-3.  (continued) 
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where  n  is  the  "chemical  potential"  and  An  is  the  change  in  chemical  species  concentration, 
so  that  the  term  p.  An  corresponds  to  the  change  in  internal  energy  due  to  dissociation.  At 
constant  volume,  we  can  simply  write  the  change  in  internal  energy  per  unit  mass  as 

AE  =  |cv,effdT  =  Cv,effAT  ,  (2.9) 

by  absorbing  the  chemical  energy  and  entropy  in  the  effective  specific  heat,  Cv,eff- 

Note  that  for  a  monatomic  gas  of  molecular  weight  M,  the  specific  heat  at  constant 
volume  per  unit  mass  is 

Cymonatomic  =  (3/2)  (R/M)  ,  (2.10a) 


where  R  (=  8.31  Joule/mole-K)  is  the  gas  constant.  For  a  diatomic  gas  (with  the  rotation 
excited  classically  and  the  vibration  not  excited)  we  have 

Cydiatomic  =  (5/2)  (R/M)  ,  (2.10b) 

while  for  a  polyatomic  gas  or  a  dissociating  diatomic  gas, 

CyPolyatomic  =  q  (R/M)  ,  (2. 10c) 

where  q  is  appropriately  large:  q  =  7/2  for  a  triatomic  molecule  with  the  vibration  not 
excited  or  for  a  diatomic  gas  with  the  vibration  fully  excited. 


For  a  chemically  reacting  gas  such  as  air  at  elevated  temperatures,  the  degree  of 
dissociation  varies  with  temperature  (and  pressure),  and  thus  the  specific  heat  ratio 

Y  =  Cp/Cv  ^2.11) 


is  no  longer  a  constant,  and  thus  Gilmore  (1955)  introduced  an  "effective  specific  heat  ratio 
for  shock  waves" 


Y’=  1+PV/E 


(2.12) 


whose  value  reverts  to  y  in  the  limit  of  a  non-reacting  ideal  gas.  For  air  at  elevated 
temperatures,  the  degree  of  dissociation  and  thus  the  values  of  y  ’  depend  on  density  as 
well  as  on  temperature. 

Figure  2-4  shows  values  of  y  ’  for  various  values  of  p/po.  Note  that  in  the  region 
where  the  O2  and  N2  molecules  dissociate,  Cv,eff  is  particularly  large  so  that  y '  is  quite 
small;  y '  =  1.1  to  1.2  corresponds  to  Cy^eff  =  (10.0  to  5.0)  (R/M)  as  against  the  much 
smaller  values  for  non-reacting  mono-  or  diatomic  gases. 
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The  structure  shown  in  Fig.  2-4  arises  because  the  various  degrees  of  freedom  are  excited 
at  increasing  temperatures  and  dissociation  and  ionization  occur  first  at  lower  densities. 
Thus,  starting  at  the  lowest  temperatures,  the  initial  decrease  in  y '  is  due  partly  to  excitation 
of  N2  and  O2  vibrational  modes,  but  even  more  to  O2  dissociation,  which  causes  large 
increases  in  internal  energy  as  the  temperature  is  increased  only  moderately.  This 
dissociation  occurs  more  abruptly  and  at  a  lower  temperature  when  the  air  density  is  lower. 
After  the  O2  dissociation  is  essentially  completed,  these  curves  level  off  or  even  rise 
slightly,  in  the  region  where  PV  increases  relatively  faster  with  increasing  temperature  than 
E  does.  At  slightly  higher  temperatures,  a  second  set  of  dips  occurs  due  to  N2 
dissociation.  Since  the  N2  dissociation  energy  is  about  twice  that  of  O2,  these  dips  occur  at 
about  twice  the  temperature  of  the  O2  dips.  At  still  higher  temperatures,  further  dips  occur 
due  to  first  ionization,  second  ionization,  etc. 

2.5  COLLISIONS  AND  NON-EQUILIBRIUM  PROCESSES 

When  energy  is  suddenly  deposited  in  a  gas,  it  is  usually  absorbed  in  a  non¬ 
equilibrium  fashion,  and  it  takes  a  finite  amount  of  time  for  this  energy  to  be  redistributed 
among  the  various  molecular  modes,  to  approach  thermal  equilibrium.  For  example,  when 
a  shock  passes  through  air,  the  energy  goes  initially  into  translational  energy  of  the 
molecules,  and  later  is  redistributed  successively  to  the  rotational,  vibrational,  electronic 
and  chemical  degrees  of  freedom.  When  air  is  bombarded  by  a  pulse  of  X-rays  from  a 
nuclear  burst,  or  by  a  pulse  of  light  from  a  laser,  the  energy  initially  goes  primarily  into 
electronic  excitation  and  ionization,  and  later  into  the  other  modes. 

This  energy  redistribution  is  carried  out  by  collisions  between  the  molecules  (and 
atoms  and  electrons,  if  present  in  significant  amounts),  and  also  perhaps  by  radiation.  The 
molecular  collision  rate  increases  with  the  gas  density.  Consequently,  near  sea  level  the 
duration  of  the  nonequilibrium  period  is  often  so  short  that  it  can  be  neglected,  and 
equilibrium  assumed  in  most  (though  not  ail)  situations.  In  contrast,  at  high  altitudes 
equilibrium  tends  to  be  the  exception  rather  than  the  rule.  However,  even  at  low  altitudes 
non-equilibrium  can  persist  if  the  energy  deposition  process  continues,  as  in  the  case  of 
delayed  beta  irradiation  from  a  nuclear  burst. 
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To  make  quantitative  estimates  of  the  equilibration  rate,  it  is  necessary  to  look  at 
details  of  the  collisional  interaction  processes.  The  gas-kinetic  collision  cross  section  of  an 
air  molecule,  a(gk)  is  in  the  range  of  1  to  3  x  1Q~15  cm^.  It  is  customary  to  define  a  gas- 
kinetic  mean  free  path  l(gk)  which  is  related  to  the  particle  number  density  n  and  the  gas- 
kinetic  collision  cross  section  o(gk)  by  the  relational 

l(gk)=l/[V2  7ina(gk)]  (2.13) 

(cf.  Jeans,  1982,  pp.  135,  145,  159). 

For  room  temperature  air  at  sea  level, n©  ~  2.55  x  lO^S  m~3  and 
l(gk)  "  6.63  X  10~8  m,  which  gives  o(gk)  ~  1.33  x  10“15  cna2 

Table  2-4  gives  some  representative  orders  of  magnimde  for  the  gas-kinetic  mean 
free  path  and  also  the  time  between  collisions  in  air  as  a  function  of  density  and  thus  of 
effective  altitude  in  the  atmosphere.  Considering  that  reaction  may  occur  on  the  average 
once  in  10^  to  10^  collisions,  we  see  that  while  equilibrium  will  tend  to  be  established 
relatively  rapidly  (fraction  of  a  second)  at  densities  greater  than  10“3  x  ambient, 
corresponding  to  altitudes  below  50  km,  yet  at  densities  below  10~^  x  ambient ,  which 
correspond  to  altitudes  above  100  km,  non-equilibrium  effects  will  predominate  to 
relatively  long  times  (tens  of  seconds). 

The  upper  atmosphere  (>  100  km)  tends  to  be  in  Photochemical  Steady  State^'^ 
rather  than  in  Local  Thermodynamic  Equilibrium  (LTE)  which  is  discussed  here. 
However,  at  the  higher  densities  or  pressures  (lower  altitudes)  that  apply  to  missile  reentry 
peak  heating  (15  to  45  km)  or  low-altitude  nuclear  fireballs  (say,  below  40  km  where  the 
pressure  is  2.9  mbar  versus  1000  mbar  at  the  surface),  LTE  is  often  a  good  first 
approximation,  although  it  fails  for  slow  processes  involving  species  present  in  very  small 
concentration. 


^  ^  Note  that  there  are  somewhat  different  definitions  for  the  mean  free  path  and  also  for  the  molecular 
collision  diameter  which  give  slightly  different  numerical  values. 

Sec,  e.g.,  U.S.  Standard  Atmosphere,  1962,  reproduced  in  the  listing  of  Physical  Constants  in 
Appendix  D. 

Photochemical  Steady  State  is  defined  as  follows.  Solar  (or  other)  ultraviolet  radiation  incident  on  the 
upper  atmosphere  produces  significant  dissociation  and  ionization,  but  does  not  raise  the  kinetic 
temperature  to  the  LTE  value  associated  with  the  ambient  dissociation  and  ionization.  See  Section  4.4 
for  a  discussion. 
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Table  2-4.  Mean  Free  Path  and  Collision  Times  as  Function  of  Density 

*  We  begin  with  Sea  Level  Atmosphere  (USAF.1965,  p.  2-19)  where  T  =  288  K  and 
Pr  =  1 .225  kg/m^  and  go  down  in  density  at  constant  temperature: 


Density  ratio 
P'Pr 

Particle  Number 
Density,  m-3 

Mean  Free  Path 
l(gk)  ,m 

Time  Between 
Collisions,  s 

1 

2.55  X  1025 

6.63  X  10-8 

1.45  X  10-1° 

10-3  (b) 

2.55  X  1022 

6.63  X  10-8 

1.45  X  10-7 

10-®(C) 

2.55x1013 

6.63  X  10-2 

1.45  X  10-^ 

10-3  (d) 

2.55  X  101® 

66.3 

0.145 

Notes: 

(a)  1 /collision  frequency. 

(b)  49  km  altitude. 

(c)  96  km  attitude. 

(d)  1 55  km  to  1 65  km,  depending  on  solar  activity-this  has  to  be  corrected  for 
temperature  530  K  to  950  K). 


For  a  chemical  reaction,  say 

A-hB<=>C-i-D  ,  (2.14) 

typically  not  every  gas-kinetic  collision  between  the  reactants  A  -i-  B  leads  to  the  products, 
C  -f-  D,  so  that  a(react)  <  o(gk).  In  a  "near-LTE"  situation  in  which  it  is  possible  to  define 
a  temperature  T,  at  least  approximately,  it  is  appropriate  to  write  down  a  reaction  rate 
coefficient  for  reaction  (2. 14)  which  is  normally  expressed  in  the  form 

k(react)  =  P  v(gk)  a(gk)  exp  { -AE/kT}  ,  (2. 1 5) 

where  v(gk)  is  a  thermal  mean  molecular  velocity  [--  10^  cm/sec,  so  that  v(gk)  a(gk) 
-  10“^0  cm^/sec].  P  is  a  probability  factor,  frequently  much  less  than  one,  while  the 
exponential  factor  represents  an  energy  barrier  AE  that  has  to  be  overcome  before  reaction 
occurs.  The  temperature  dependence  of  the  reaction  rate  coefficient  comes  mainly  from  the 
exponential  term. 


For  the  reaction  (2.14),  where  species  concentrations  are  nx  (mol./cm^; 
X  =  A,  B,  C,  D)  and  the  forward  and  backward  rate  coefficients  are  respectively  kf<T)  and 
kb(T),  in  a  steady  state  the  number  of  forward  reactions  per  unit  volume  and  unit  time, 
kf(T)  nA  nfi,  is  equal  to  the  number  of  back  reactions  per  unit  volume  and  unit  time,  kb(T) 
nc  nD-  This  is  the  Principle  of  Detailed  Balancing,  which  states  that  when  equilibrium  is 
reached  in  a  reaction  system,  any  particular  chemical  reaction  and  its  exact  reverse  must,  on 
the  average,  be  occurring  at  the  same  rate.  It  gives  a  useful  relation  between  the 
equilibrium  constant  Keq(T)  for  reaction  (2.14)  and  the  rate  constants  kf(T)  and  kb(T) 

Keq(T)  =  nc  no/nA  nB  =  kf<T)/kb(T)  .  (2.16) 

An  interesting  and  important  example  of  a  non-equilibrium  process  in  high- 
temperature  air  is  provided  by  the  Nitric  Oxide  Freeze-In  (see  footnote  4  on  p.  2-4). 
Reference  to  Fig.  2-3  or  -  standard  sources  shows  that  at  temperatures  above  2000  K, 
equilibrium  air  cont?‘n'  '.gnificant  quantities  of  NO  (1  to  10  mole-percent,  or  half  that  in 
units  of  particles  per  air  atom).  The  NO  is  produced  by  a  variety  of  chemical  reactions, 
with  an  important  contributor  being  O  -t-  N2  +  M  (=  some  third  body);  it  disappears  largely 
by  the  reaction 

NO NO  *>  N2  +  O2  ,  (2.17) 

which  is  exothermic  by  1.9  eV,  but  has  a  large  activation  energy  AE  =  5.9  eV 
(135  Kcal/mole),  and  thus  becomes  exceedingly  slow  at  temperatures  below  2000 
Accordingly,  when  high-temperature  air  cools,  the  NO  remains  for  a  long  time  in  above¬ 
equilibrium  concentrations.  It  forms  NO2  which  is  an  irritating  oxidant  itself,  and  which 
also  reacts  photochemically  with  ambient  atmospheric  hydrocarbons  to  form  ozone  and 
other  oxidants,  so  that  this  process  is  an  important  source  of  photochemical  smog  arising 
from  the  products  of  combustion. This  freeze-in  process  is  important  in  the  combustion 
products  from  heat  engines,  in  the  wakes  behind  reentering  objects,  and  in  nuclear 
fireballs. 

If  the  chemical  composition  in  a  given  situation  cannot  be  described  by  near¬ 
equilibrium,  in  terms  of  an  effective  temperature,  the  following  points  should  be  noted: 


Because  the  rate  constant  -  exp  -(AE/kT)  =  exp  -(5.9/0. 18)  =  e“^^  ~  5  x  10"^^. 

'  ^  In  the  stratosphere,  due  to  the  lower  hydrocarbon  concentration  and  greater  ultraviolet  intensity  than 
near  the  ground,  NO  usually  destroys  rather  than  {xoduces  ozone. 
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1 .  One  improvement  which  is  ftequently  useful  is  to  describe  different  degrees  of 
freedom  (e  g.,  molecular  vibration,  dissociation,  or  electronic  excitation  )  in 
terms  of  their  own  effective  temperature(s)  Tjit  which  may  not  always  be  equal 
to  the  gas-kinetic  (i.e.,  translational)  temperature. 


2. 


3. 


For  high-temperature  rate  constants,  see  the  evaluated  data  collections  of 
Baulch  et  al.,  1972,  1973,  1976,  for  near-ambient  temperatures,  Baulch  et  al., 
1984,  and  for  charged  species,  DNA,  1972. 

Sets  of  coupled  chemical  reactions  are  numerically  frequently  difficult  to 
analyze  because  as  the  system  approaches  a  steady  state,  small  fluctuations  in 
concentration  about  the  steady-state  value  tend  to  be  amplified  by  numerical 
techniques.  This  "stiff"  behavior  can  be  avoided  by  using  a  procedure 
originally  due  to  Gear,  1971. 


2.6  SOME  COMMENTS  ON  HEAT  TRANSFER  IN  HIGH- 
TEMPERATURE  AIR 

In  the  present  high-temperature  regime,  heat  transfer  is  a  particularly  complex 
subject.  For  a  discussion,  see,  e.g..  Bond  et  al.,  1965,  Penner  and  Olfe,  1968,  and 
Zeldovich  and  Raizer,  1966,  1967,  and  on  the  aerodynamics,  see  Anderson,  1982,  1984, 
Gazley,  1961,  and  Martin,  1966. 

Heat  transfer  occurs  by  conduction  (molecular  transport),  by  convection  and  by 
radiation.  Convective  heat  transfer  is  normally  faster  than  conductive  transfer  in  gases,  and 
predominates  at  temperatures  below  several  thousand  K.  At  higher  temperatures  radiative 
heat  transfer  is  much  more  rapid.  Typically,  in  reentry,  convection  dominates  radiation  up 
to  several  thousand  degrees  K;  thus,  for  ordinary  earth  satellite  or  ICBM  entry  where  the 
entry  velocity  is  t>'pically  7  km/sec  and  the  maximum  temperature  is  ~  7500  K,  the  effects 
of  radiative  heat  transfer  are  normally  not  very  important,  but  at  temperatures  of  12000  K 
or  higher,  which  correspond  to  entry  velocities  of  perhaps  11  km/sec,  radiative  heat 
transfer  dominates.  In  nuclear  fireball  problems  also  convection  dominates  radiation  up  to 
a  few  thousand  degrees. 

Note  that  : 

1 .  In  the  temperature  regime  considered  here,  a  large  fraction  of  the  thermal 
energy  is  transported  as  heat  of  dissociation  of  the  O2  and  N2  molecules  rather 
than  as  kinetic  and  rotational  energy  of  the  molecules. 

2.  Convective  heat  transfer  increases  very  much  when  the  flow  goes  turbulent. 
The  determination  of  when  a  flow  goes  turbulent  and  how  large  the  heat 
transfer  then  becomes  is,  in  general,  complex  and  very  much  dependent  on  the 
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detailed  character  of  the  flow,  so  that  a  simple  and  useful  account  of  this 
subject  cannot  be  presented  here.i^ 

3 .  Air  does  not  radiate  at  all  like  a  black  body  because  there  are  many  absorption 
bands  interspersed  by  transmissive  regions.  Under  reentry  conditions  (T  ~ 
5000-7000  K)  radiative  heat  transfer  sometimes  increases  very  much  more 
rapidly  with  temperature  than  blackbody  radiation;  S.S.  Penner  (private 
communication)  mentions  that  under  certain  conditions  the  radiative  heat 
transfer  is  proportional  to  as  against  for  blackbody  radiation. 


'  ^  See  Ch.  5  in  Part  II  of  these  Notes. 
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3.  RADIATION  FROM  HIGH-TEMPERATURE  AIR 


3.1  INTRODUCTION 

The  radiation  from  air  and  from  other  gases  at  high  temperatures  is  highly  specific 
spectrally,  and  one  frequently  needs  some  very  detailed  information  about  the  particular 
chemical  species  and  radiating  systems  that  are  important  At  one  atmosphere  pressure  and 
below,  atomic  spectral  lines  of  the  atmospheric  gases  tend  to  have  a  very  narrow  spectral 
line  width  (<  1  cm~l),  so  that  a  given  air  mass  is  likely  to  be  optically  thick  at  the  centers  of 
the  stronger  lines  and  optically  thin  in  other  spectral  regions.  This  gives  rise  to  deviations 
from  Beer's  Law,  in  that  the  attenuation  at  wavelength  X  in  path  length  x  cannot  always  be 
written  as  exp(-kxx)  where  kx  depends  on  X  but  not  on  x.^"^  Molecular  transitions  involve 
bands  made  up  of  hundreds  or  thousands  of  spectral  lines,  but  conceptually  they  are  not  too 
difficult  to  describe.  In  a  diatomic  molecule,  either  one,  two,  or  all  three  of  the  quantum 
parameters  designating  the  rotational,  vibrational,  and  electronic  state  of  the  molecule  can 
change.  These  changes  typically  cause  emission  or  absorption  in  the  microwave,  infrared, 
and  visible-ultraviolet  spectral  regions,  respectively.** 

For  orientation,  note  that  a  1  cm  thick  layer  of  air  at  4000  K  and  1  atmosphere 
pressure,  such  as  in  a  reentry  vehicle  boundary  layer,  is  quite  thin  optically  and  thus 
radiates  very  much  less*^  than  a  blackbody  at  this  temperature,  but  a  1000-m  layer  of  air  at 
this  temperature  and  pressure,  as  in  a  nuclear  fireball,  is  optically  thick  and  radiates  like  a 
blackbody  .2® 


*  ^  This  issue  is  discussed  in  more  detail  in  Part  II,  Section  7.2.3  in  the  context  of  the  LOWTRAN  code. 

*  *  Triatomic  molecules  behave  similarly,  except  that  (if  linear)  they  have  one  (doubly-degenerate) 

rotational  mode  and  three  vibrational  modes  (one  doubly-degenerate),  while  if  nonlinear  they  have  three 
rotational  and  three  vibrational  modes,  all  nondegenerate. 

*  ^  By  a  factor  10~*-see,  e.g..  Fig.  3-10. 

At  temperatures  in  the  1000  K  to  10,000  K  range  both  a  near-gray  body  and  an  air  parcel  emits  most  of 
its  radiation  in  the  "near- visible"  spectral  range,  say  0.2  pun  to  2  pm.  This  is  so  because  many 
absorption  bands  are  located  in  this  spectral  range  in  which  most  of  the  energy  of  a  Planck  (blackbody) 
spectrum  lies  in  the  temperature  range  under  consideration.  In  general,  this  spectral  range  or  the  visible 
is  considered  here. 


We  begin  in  Section  3.2  by  discussing  blackbody  radiation  (which  gives  the 
maximum  radiance  at  a  given  temperamre)  from  the  standpoint  of  applications  to  the 
radiation  from  solid  surfaces.  Then,  in  Section  3.3,  we  discuss  the  equation  of  radiative 
transfer,  including  effects  of  scattering  as  well  as  of  emission  and  absorption.  The  limits  of 
optically  thick  and  optically  thin  media  are  discussed,  and  Rosseland  and  Planck  means 
(which  are  appropriate  for  these  two  limits)  are  introduced.  Section  3.4  treats  atomic  line 
radiation  to  point  out  some  generic  characteristics  of  absorption  cross  sections,  line  width 
and  line  shapes,  and  saturation  in  terms  of  "Curves  of  Growth."  Section  3.5  analyzes 
molecular  band  radiation  for  the  particular  example  of  high-temperature  air  where  six 
molecular  band  transitions  (O2  Schumann-Runge,  NO  [3  and  y,  N2  First  Positive  and 
Second  Positive  and  N2'^  First  Negative)  have  frequently  been  used  to  describe  the  bulk  of 
the  near-visible  thermal  radiation  which  predominates  in  the  temperature  range  of  concern 
here.  Section  3.6  discusses  the  continuum  ("free-free"  and  "free-bound")  radiation  which 
contributes  especially  at  the  higher  temperatures. 

Section  3.7  uses  these  results  to  present  the  absorption  coefficient  of  high- 
temperature  equilibrium  air.  Section  3.8  applies  these  results,  indicating  in  quantitative 
terms  that  in  the  present  temperature  range  (1,(X)0  K  to  10,000  K)  radiation  from  the  air  is 
much  less  than  blackbody  radiation  at  the  same  temperature. 

It  must  be  stressed  that  the  present  discussion  is  not  a  substitute  for  a  course  in 
atomic  and  molecular  spectroscopy  and  quantum  mechanics,  which  would  explain  how  the 
electronic  structme  is  built  up^^  and  how  one  calculates  the  transition  probabilities  of  the 
various  transitions.  For  this,  see  for  instance  Condon  and  Shortley,  1935,  Goody,  1964, 
Herzberg,  1937,  1944,  1945,  1950,  1966,  etc.,  Penner,  1959,  Slater,  1960,  1963,  and 
Zeldovich  and  Raizer,  1966,  1967.  Here  we  try  to  sketch  out  the  underlying  structme  in 
qualitative  terms  to  explain  how  the  various  pieces  of  the  problem  fit  together,  present  some 
key  analyses,  and  show  numerical  values  for  some  key  optical  transitions  to  give  the  reader 
a  feeling  for  the  orders  of  magnitude. 

3.2  BLACKBODY  RADIATION 

A  blackbody  is  defined  as  a  medium  that  absorbs  all  the  incident  radiation  that  falls 
on  it.  Since,  in  equilibrium,  the  amount  of  energy  absorbed  must  equal  the  amount 
emitted,  the  equilibrium  radiant  energy  emitted  from  unit  area  of  a  blackbody  in  unit  tim 

^ '  Appendix  A  gives  a  very  brief  discussion  of  spectroscopic  notation,  which  is  used  extensi  ■'cly  here. 


represents  an  upper  limit  for  the  emission  of  radiant  energy  from  any  substance  which  is  at 
the  same  temperature  as  the  blackbody.  The  radiation  from  a  blackbody  at  a  given 
wavelength  X  or  frequency  v  depends  only  on  X  (or  v)  and  on  T,  and  gives  the  maximum 
radiance  from  any  source  at  this  temperature.  The  blackbody  spectral  radiance  B(A,,T)  is 
given  by  the  following  expression: 

B(X,T)  =  Cl  >.-5  [exp(c2AT)-l]-l  (W/m2-sr-^im)  ,  (3.1) 

where  (for  X  in  |im,  T  in  K): 

Cl  =  2hc2  =  1.191  X  10^  (W-|i.m5/m2-sr-|im) 

C2  =  hc/k  =  14390  (K-|im)  . 

Figures  3-1  and  3-2  show  the  spectral  radiance  at  different  temperatures  as 
functions  of  wavelength  X  (pm)  and  wavenumber  v  (cm~^).  Note  that  the  symbol  v  is 
sometimes  used  for  frequency  (Hz)  and  sometimes  for  wavenumber  (cm“^).  in  the  present 
document  we  normally  use  wavelength  rather  than  frequency  or  wavenumber  as  spectral 
variable,  except  in  some  particular  discussions  such  as  that  in  Section  3.3.3  and  in 
Appendix  C.  Overall,  the  notation  used  in  the  literature  is  very  confusing.  For  instance, 
many  people  to  refer  to  W/cm^  rather  than  to  W/m2.  See,  e.g.,  IR  Handbook,  1985,  Ch. 
2.  A  good,  concise  summary  is  given  in  USAF,  1965,  Appendix  B,  which  is  the  source  of 
Figs.  3-1  and  3-2  . 

If  one  integrates  the  spectral  radiance,  B(X,T),  one  obtains  the  total  radiant  energy 

density 

u  =  a  T^  (a  =  87c5k4/15c3h3)  (J/m^)  (3.2a) 


or  the  radiant  exitance  of  a  blackbody 
E  =  (c/4;t)aT4  =  a  t4  , 


(W/m2) 


(3.2b) 


where  the  Stefan-Boltzmann  constant  a  =  5.670  x  10"^  W/m^K^. 


Differentiating  Eq.  (3.1)  with  respect  to  wavelength  and  equating  3B(X,T)/3X  to 
zero  gives  Wien's  displacement  law,  which  states  the  relation  between  Xmax.  the 
wavelength  at  which  the  spectral  radiance  B(X,T)  is  a  maximum,  and  temperature 

Xfnax  T  =  2.898  X  10^  (pm  K)  .  (3.3) 
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RADIANCE  per  WAVENUMBER  I W  m'*' sr“' (cm' 


Figure  3-2.  Spectral  Radiance  of  a  Blackbody  as  Function  of  Wavenumber. 
•  (Source:  USAF,  1965,  Appendix  B) 
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Actual  solid  surfaces  often  radiate  proportionally  to  a  blackbody  with  an  effective 
emittance22  e  =  e(X),  so  that  the  spectral  radiance  L(X,T)  (W/m2-sr-|im)  from  a  solid 
"Lambertian"  surface^^  is  given  by  the  expressions^ 

L(X,T)  =  £a,T)B(>.,T)  .  (3.4) 

Figure  3-3  gives  the  reflectance 

reflectance  =  1  -  emittance  (3.5a) 

of  various  natural  (solid  or  liquid)  materials  at  room  temperature  as  a  function  of 
wavelength.  If  there  is  a  general  rule,  it  is  that  many  condensed  phase  materials  have  a 
large  emittance  (small  reflectance)  at  wavelengths  above  5  jim  to  10  |im,  i.e.,  at  frequencies 
below  the  Debye  frequency  of  lattice  vibrations.  In  the  visible,  transparent  materials  have  a 
very  low  emittance,  but  in  the  UV  the  emittance  sometimes  increases  again  due  to  atomic 
transitions.  The  temperature  dependence  of  the  emittance  of  most  materials  is  not  strong, 
except  when  there  is  a  phase  change  or  other  structural  transition. 

An  important  principle  is  Kirchoff  s  Law:  at  each  wavelength 

emittance  =  absorptance=  1  -  reflectance  (3.5b) 

which  is  universally  applicable. 


If  e  =  constant,  independent  of  X,  this  defines  a  gray  body. 

A  Lambenian  surface  is  one  which  reflects  and  emits  energy  diffusely. 

Some  authors  use  "spectral  irradiance"  E(X.T)  (W/m^-pm)  which  is  the  integral  of  L(X,T)  over  solid 
angle-see,  e.g., //?  Handbook,  1985,  Ch.l. 
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REFLECTANCE 


Figure  3*3.  Typical  Reflectance  of  Various  Natural  Surfaces. 
(Source:  McCiatchey  et  al.,  1972) 
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3.3  EQUATION  OF  RADIATIVE  TRANSFER 
3.3.1  Without  Scattering 

The  equation  of  transfer  describes  the  change  in  the  intensity  function  I  =  I(X,x) 
where  x  =  distance,  in  a  medium  characterized  by  an  absorption  coefficient  ka  (X)  (m“*): 

ka(X)  =  Ij  nj  aa(j;X)  .  (3.6) 

Here  nj  =  number  density  of  absorbers  of  species  j,  whose  absorption  cross  section  at 
wavelength  X  is  aa(j;X).  In  fact,  this  expression  for  the  absorption  coefficient  kg  has  to  be 
corrected  for  stimulated  emission  (see  Appendix  C)  to  give 

k'a(X)  =  ka(X)  [1  -  exp(-C2AT)]  (3.7) 

(recall  that  C2A-T  =  hv/kT) . 

For  the  transport  of  radiant  energy,  it  is  customary  to  write  down  the  Equation  of 
Transfer,  This  gives  the  radiant  intensity  I  =  I(X,x)  as  a  function  of  an  appropriate  source 
function  S*(X,x)  and  of  the  absorption  coefficient  k'a(X)  (see,  e.g.,  Aller,  1963,  Ch.4  and 
Ch.5)  as 

dI(X,x)/dx  +  k'a(X)  I(X,x)  =  k'a(X)  S*(X,x)  (3.8) 

and  for  a  medium  in  which  we  define  a  temperature  at  each  point,  T  =  T(x),  we  have 

S*(X,x)  =  B[X,T(x)]  ,  (3.9) 

so  that 

dI(X,x)/dx  =  k'a(X)  { B[X,T(x)]  -  I(X,x) }  .  (3. 10) 

Let  us  now  demonstrate  how  blackbody  radiation  relates  to  the  radiation  from  a 
uniform  layer  of  gas  at  temperature  T  and  of  thickness  x.  If  one  uses 

u  =  k'a  X  (3.11) 

as  an  independent  variable  and  suppresses  X,  the  transfer  equation  (3.10)  becomes 

dl/du  + 1  =  B  or  (d/du)(Ie“)  =  B  e“  .  (3.12) 

Now,  with  the  boundary  condition 

l(x  =  0)  =  0  ,  (3.13) 
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i.e.,  no  incident  radiant  from  outside  our  volume  of  interest,  we  find 


I(x)  =  B(1  -  e-'^'ax)  ,  (3.14) 

so  that  for  the  two  limits  k'ax  «  1  and  »  1  we  have 
Optically  thin,  k’ax  «  1,25 

I(x)  =  Bk'aX  ,  (3.15) 

so  that  the  radiation  is  proportional  to  number  of  atoms,  and 
Optically  thick.  k'aX  »  1,26 

I(x)  =  B(X,T)  ,  (3.16) 

i.e.,  the  gas  radiates  like  a  blackbody. 

The  next  problem  is  to  determine  ka(X),  the  absorption  coefficient,  and  this  is  done 
in  the  next  three  sections,  respectively,  for  atomic  line  radiation  (Section  3.4),  molecular 
band  radiation  (Section  3.5),  and  continua  (Section  3.6).  In  these  three  sections 
3.4  to  3.6  we  evaluate  Oa(j;X);  the  effects  of  number  density  or  of  optical  thickness  are 
taken  up  in  Section  3.7  thereafter.  However,  before  we  do  this,  let  us  write  down  the 
equation  of  transfer  in  the  presence  of  scattering  particles  in  Section  3.3.2,  and  set  this 
discussion  of  radiative  transfer  in  context  with  a  general  discussion  of  opacity  coefficients 
in  Section  3.3.3. 

3.3.2  With  Scattering 

If  one  deals  only  with  atolHs  and  molecules  at  reasonably  low  density  and  at 
relatively  long  wavelengths,  there  is  very  little  scattering  and  thus  the  equation  of  transfer  in 
the  form  (3.8)  is  satisfactory.  However,  if  there  are  particles  of  radius  a  such  that  the  Mie 
parameter22 

q  =  2TcaA  (3.17) 


25  Which  is  gcxxl  to  10%  for  k'aX  <  0.2. 

26  Which  is  good  to  10%  for  k'aX  >  1.5. 

22  In  1908  G.Mie  gave  a  rigorous  solution  for  the  diffraction  of  a  plane  electromagnetic  wave  by  a 
homogeneous  sphere--see,  e.g..  Bom  &  Wolf,  1980,  or  Van  de  Hulst,  1957. 


is  not  very  much  less  than  unity,  then  scattering  may  be  as  important  as  absorption  or  even 
more  so.  In  general,  q  «  1  defines  the  Ravleigh  scattering  re^ime.^^  q  ~  1  is  called  the 
Mie  region,  and  q  »  1  is  the  domain  of  Geometrical  Optics.  Scattering  effects  will  be 
discussed  in  some  detail  in  Chapter  7  of  Part  H,  to  be  prepared  later. 

If  there  are  n  particles  per  unit  volume,  each  of  absorption  cross  section  Oa  and  of 
scattering  cross  section  Os  =  Os(0),  then  we  can  define  an  absorption  coefficient  ka(X)  from 
Eq.  (3.6)  and  a  scattering  coefficient 

ks(0)  -  n  dCTs/de  (3.18) 

where  dOs/d0  is  the  differential  scattering  cross  section.29 

The  equation  of  transfer  describes  conservation  of  total  radiant  intensity  in  a 
volume  element,  when  account  is  taken  of  absorption  and  emission  (terms  in  k'a  and 
B(X.,T))  as  well  as  of  scattering  (terms  in  kg).  We  write  down  the  equation  of  transfer  for 
the  intensity  function  I  =  I(x,  m  (p)  in  a  semi-infinite  slab  of  uniform  material,  characterized 


by  X  >  0,  in  a  direction  (0,9):^® 

(ks+ka')l-iuil/dx  = 

(1/47c)|J  dp'd(p'ks(©)  I(x,p’,(p’)  +  ka’  B(T) 

(3.19a) 

l^s  = 

(l/4jr)//dp'  dtp'  ks(0) 

(3.19b) 

COS  0  = 

cos  0  cos  0'  +  sin  0  sin  0’  cos  (tp-tp') 

(3.19c) 

= 

cos  0  ;  p'  =  cos  0'  . 

(3.19d) 

This  equation  is  written  down  here  to  indicate  that  absorption  and  scattering  enter  in 
somewhat  different  ways.^^  By  comparing  Eq.  (3.19)  with  (3.8)  we  see  that  scattering 
and  absorption/emission  enter  in  somewhat  different  ways.  Thus  k'a-but  not  ks—is 

corrected  for  stimulated  emission.  The  difference  in  sign  of  the  term  dl/dx  depends  on  the 


Lord  Rayleigh  explained  the  blue  of  the  sky  as  due  to  the  scattering  of  sunlight  by  air  molecules. 
Molecules  are,  of  course,  extremely  small  (a  ~  3  x  10-8  cm)  so  that  q  «  1  and  thus  the  scattering  per 
molecule  is  very  small,  but  the  optical  thickness  of  a  path  through  the  atmosphere  is  large  because 
there  are  so  many  molecules  in  the  path.  If  q  «1,  the  scattering  cross  section  per  particle  is 
proportional  to  \fK^,  so  that  blue  light  is  scattered  much  more  than  red  light,  which  explains  why  the 
sky  is  blue.  (This  is  a  great  oversimplification-see,  e.g.,  Minnaert,  1954,  p.  238f,  and  Van  de  Hulst, 
1957,  p.  414.).  Also,  the  sun  is  red  at  sunrise  and  sunset  because  it  is  seen  through  a  long  path  in  the 
dusty  lower  atmosphere,  which  ineferentially  removes  the  blue  from  the  path. 

Note  that  Cj  and  k,  depend  on  wavelength  \  as  well  as  on  angle  9. 

See,  for  instance,  Chandrasekhar,  1960,  Ch.l,  or  Bond,  Watson  and  Welch,  1%5,  Ch.lO  . 

Conventionally,  p  =  cos  0  is  defined  so  that  p  >  0  corresponds  to  radiation  coming  out  of  the  material 
and  p  <  0  to  radiation  going  into  the  material . 
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sign  convention  used  for  p..  The  right-hand  side  of  Eq.  (3.19a)  contains  an  emission  term 
k’aB  and  a  scattering  term  (l/47c)//d|j.'  dtp'  ks(0)  I(x,|i',(p'),  while  the  scattering  term  kg 
has  been  integrated  over  all  scattering  angles.  We  shall  not  discuss  the  solution  of  the 
equation  of  radiative  transfer  (3.19),  which  becomes  very  complicated  (see,  e.g.,  Aller, 
1963,  Chandrasekhar,  1960,  Goody,  1964,  or  Kondratyev,  1969),  but  merely  present  it 
here  to  indicate  how  the  introduction  of  scattering  effects  complicates  the  problem. 

3.3.3  Opacities^^ 

In  an  optically  thin  ("transparent")  medium  we  ask  how  much  the  intensity  of  a 
beam  of  incident  radiation  is  reduced,  i.e.,  we  want  to  know  the  (relatively  small) 
absorption  coefficient  for  a  situation  in  which  Beers'  Law  generally  holds.  By  contrast,  in 
an  optically  thick  ("opaque")  medium  the  critical  question  is  how  far  a  photon  travels,  or 
what  is  the  mean  free  path.  This  concept  of  optically  thin  versus  optically  thick  media  is 
straightforward,  but  its  evaluation  is  not,  because  of  the  complications  arising  from  the 
various  physical  factors  in  the  equation  of  transfer. 

As  has  been  shown  above,  in  the  presence  of  scattering  as  well  as  absorbing 
particles,  the  equation  of  transfer  is  an  integro-differential  equation  which  is  normally 
solved  in  a  number  of  approximations.  It  is  customary  (see,  for  example,  Bond,  Watson 
and  Welch,  1965,  Chapter  10)  to  distinguish  between  the  following: 

a.  Diffusion  approximation,  obtained  by  setting  I  ~  B  in  the  transfer  equation 
(3.10).  This  is  valid  for  optically  thick  regions  where  photons  emitted  in  the 
region  have  a  high  probability  of  reabsorption.  Here  one  introduces  the 
Rosseland  mean  opacity  ,  and  because  we  are  talking  about  an  optically  thick 
region,  it  is  useful  and  customary  to  use  the  Rosseland  mean  free  path. 

b.  Emission  approximation,  obtained  by  setting  I  «  B  in  the  transfer  equation 
(3.10).  This  is  valid  for  optically  thin  regions  from  which  photons  emitted  in 
the  region  can  escape  freely,  so  that  scattering  is  not  important.  Here  one 
introduces  the  Planck  mean  opacity. 

c.  In  the  absence  of  emitters,  there  is  a  cold  non-radiating  medium  through  which 
radiation  simply  diffuses.  We  shall  not  discuss  this  case  further  here,  but 
Chapter  7  in  Part  II  will  treat  various  applied  problems  relevant  to  the  earth's 
atmosphere. 


The  tenn  "opacity"  is  used  technically  for  the  absorption  coefficient  integrated  over  the  whole  spectrum 
with  a  blackbody  weighting,  cither  B(v,T)  or  B(X,T). 


1 .  Diffusion  Approximation  in  an  Optically  Thick  Medium 

In  an  optically  thick  region  one  introduces  an  absorption  mean  free  path^^ 

la(v)  =  l/ka(v)  .  (3.20) 

where  the  absorption  coefficient  ka(v)  is  defined  by  Eq.  (3.6),  and  a  total  mean  free  path^ 

l/l’(v)  =  ka(v)  +  ks(v)  .  (3.21) 

where  the  correction  for  stimulated  emission  only  occurs  in  the  absorption,  not  the 
scattering  coefficient.  To  introduce  a  frequency-averaged  form  of  the  transfer  equation 
(3.10)  or  (3.19),  we  introduce  the  Rosseland  mean  free  path  1r  (or  Rosseland  mean 
absorption  coefficient  "opacity"  kR)  as 

Ir  =  l/kR  ={/  l'(v)  aB(v,T)/aT  dv}  /  {J  aB(v,T)/aT  dv}  .  (3.22) 

The  physical  meaning  of  this  is  the  following.  In  an  atmosphere  with  a  spatial  gradient 
dl(v)/dx,  we  define  a  Rosseland  mean  opacity  as^^ 

1/kR  =  {Ir(v)  [dl(v)/dx]  dv}  /  {/[dl(v)/dx]  dv}  ,  (3.23) 

replace  I(v)  by  B(v)  ( hypothesis  for  case  1)  and  put 

dB(v,T)/dx  =  [dB(v,T)/dT].[dT/dx]  ,  (3.24) 

which  leads  to  the  result  (3.22). 

2 .  Emission  Approximation  in  an  Optically  Thin  Medium 

In  an  optically  thin  medium,  the  temperature  will  often  be  uniform  so  that  dVdx  or 
dB/dT  is  not  an  appropriate  quantity,  and  indeed  scattering  is  unimportant  (because  the 
medium  is  optically  thin).  In  this  case  I(v)  «  B(v),  and  one  asks  for  the  total  flux  from 
the  surface,  Jp,  which  in  this  ("Planck")  limit  is  given  by  the  multiple  integral 

Jp  =  lo**  la  1 0  B(v,T)  k’a(v)  dv  dQ  dx  ,  (3.25) 


As  is  conventional,  in  the  present  theoretical  discussion  we  describe  the  spectral  dependeiKe  in  terms  of 
frequency  v  rather  than  wavelength  X. 

The  scattering  coefficient  kg  used  here  is  the  total  scattering  coefficient,  integrated  over  all  scattering 
angles  6.  (The  confusing  notation  is  conventional!). 

The  Rosseland  mean  is  frequently  used,  but  "it  has  little  obvious  physical  significance  throughout  the 
upper  visible  layers  of  stellar  atmosphere"  (Aller,  1963,  p.  239). 
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where  dx  =  ds  cos0  is  distance  normal  to  the  surface  of  the  slab^^  and  the  integration  over 
I  the  solid  angle  Q  is  over  outward  directions  only.  In  an  optically  thin  medium  the 

temperature  is  essentially  constant,  so  that  if  one  evaluates  the  spatial  integral  over  dQ  and 
dx  in  Eq.  (3.25),  one  obtains  the  result 

Jp  =  (27ta)  Ax  Jka'(v)  B(v)  dv  =  27C  a  kp  Ax  a  ,  (3.26) 

where  "a"  is  a  numerical  factor^”^  and  the  Planck  mean  opacity  kp  is  given  by^* 

kp  =  /  k'a(v)  B(v,T)  dv  /  J  B(v,T)  dv  .  (3.27) 

Note  that  the  expression  (3.27)  for  kp  looks  rather  different  from  the  expression 
f  (3.23)  for  kp,  but  if  we  express  the  Rosseland  mean  in  terms  of  J(v),  the  actual  net  flux,^^ 

then  since 

J(v)  =  [l/k’(v)l  dB(v,T)/dT  ,  (3.28) 

^  we  have,  in  analogy  to  Eq.  (3.26), 

1/1r  =  kp  =  /  k’(v)  J(v)  dv  /  J  J(v)  dv  .  (3.29) 

The  difference  in  weighting  of  Planck  and  Rosseland  means  can  give  significantly  different 
numerical  values  at  a  given  temperature  and  density,  and-as  has  been  indicated  at  the 
^  beginning  of  this  Section--they  are  used  under  different  conditions. 

Table  3-1  lists  numerical  values  for  kp  (cm~l)  and  Ip  (cm)  for  high- temperature  dry 
air,  for  a  range  of  temperatures  and  densities  of  interest  in  the  present  application.  Note 
that  frequently  kp  and  1/lp  do  not  vary  monotonically  with  the  density  ratio  p/po  and  the 
temperature.  While  this  comes  in  part  from  the  fact  that  different  species,  of  different 
absorption  coefficients,  are  present  at  different  density  ratios  at  the  same  temperature,  yet  in 
part  this  non-linear  behavior  comes  from  the  saturation  of  the  curves  of  growth  of  the 
individual  spectral  lines  which  is  discussed  in  Section  3.4  below— see  in  particular  Fig.  3-5. 


^  ^  s  =  distance  along  the  path  of  the  radiation  and  0  =  angle  relative  to  normal. 

37  Note  that  the  slab  is  not  optically  thin  for  angles  9  such  that  ka(v)  Ax  /  cos9  >  1,  and  thus  the 
numerical  factor  a  in  Eq.  (3.26)  is  equal  to  1 .8  rather  than  2,  as  the  direct  analysis  would  suggest. 
Recall  that  J  B(v.t)  dv  =  a  T*. 

In  an  optically  thick  medium  the  temperature  varies,  so  that  J(v)  depends  on  the  temperature  profile 
within  the  medium. 
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Table  3*1.  Planck  Mean  Absorption  Coefficient  and  Rosseland  Mean  Free  Path 
for  Air  at  Different  Temperatures  and  Densities 
(Source:  Landshoff  and  Magee,  1969) 


A.  Planck  Mean  Absorption  Coefficient,  kp(cm~'') 


p/po 

1 

10-3 

10-6 

T=  2000  K 

1.17E-4 

1.15E-7 

1.01E-10 

4000  K 

8.04E-4 

4.70E-8 

3.27E-12 

6000  K 

5.79E-3 

4.83E-7 

1.1  IE-11 

9000  K 

3.51  E-5 

1.13E-5 

3.74E-9 

12000  K 

5.25E-1 

3.37E-4 

2.72E-9 

B.  Rosseland  Mean  Free  Path,  Ir  (cm) 


p/po 

1 

10-3 

10-6 

2000  K 

2.06E+12 

2.06E+15 

2.05E+18 

4000  K 

1.45E+4 

2.75E+8 

3.63E+12 

6000  K 

6.41  E+2 

1.03E+7 

5.00E-h11 

9000  K 

5.41  E+1 

9.07E+5 

3.98E+9 

12000  K 

1.29E+1 

5.14E+4 

5.63E+9 
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3.4  ATOMIC  LINE  RADIATION:  TRANSITIONS,  f-NUMBER 
AND  ABSORPTION  CROSS  SECTION,  LINE  SHAPE, 

CURVES  OF  GROWTH 

Line  radiation  provides  a  good  introduction  to  a  number  of  physical  concepts  that 
are  very  important  in  their  application  to  the  radiation  from  high-temperature  air. 

As  a  simple  example,  consider  sodium  vapor  at  low  pressure.  As  the  temperature 
increases,  the  vapor  looks  yellow  because  it  radiates  strongly  in  the  Na  doublet  589  nm  to 
589.6  nm.  This  is  the  transition:  (cf.,  Radzig  and  Smirnov,  1985,  p.  225) 


• 

X  (nm) 

Aui(108  s-1'  ) 

fiu 

gu/gl 

32s  1/2 

-  32p,n 

589.6 

0.61 

0.32 

2/2 

■  32P3/2 

589.0 

0.61 

0.64 

4/2 

Here  X.  is  the  wavelength  of  the  transition,  Aui  is  the  Einstein  transition  probability 
("A-  coefficient")  for  the  transition,  fiu  is  the  "f-number"  or  oscillator  strength,  and  gjg\  is 
the  ratio  of  statistical  weights  of  upper  and  lower  states  for  this  multiplet  transition.  These 
quantities  will  be  defined  and  explained  below.  By  adding  ih&f -numbers,  fiu,  of  these  two 
transitions  we  find  f  =  0.96  =  1 

Electronic  transitions  involve  changes  in  quantum  numbers,  and  there  are  selection 
rules  which  determine  the  relative  probability  of  a  given  transition.  In  atoms,  the  most 
probable  transitions  are  those  for  which  the  total  orbital  angular  momentum  L  changes  by 
one  unit,  i.e.,  AL  =  ±  1,  while  the  spin  angular  momentum  S  does  not  change.  These  are 
called  dipole-allowed  transitions  and  arc  traditionally  characterized  by  an  oscillator  strength 
or  f-number,'^^  fiu.  For  such  a  dipole- allowed  transition,  first-order  permrbation  theory 
gives  the  result 

fiu  =  (167i3mc/3e2h)  v(Hz)  l(ulrll)|2  ,  (3.30a) 

where  the  quantity  (ulrll)  is  the  dipole  transition  matrix  element  for  the  electron  coordinate  r. 
In  terms  of  the  more  general  Einstein  A-coefficient,  we  have 

fiu  =  1.349  X  1021  Aui/v(u1)2  .  (3.30b) 


This  means  that  in  terms  of  the  Old  Quantum  Theory  (Bohr,  1914  pre-Quantum  Mechanics)  the 
Na  D-line  corresponds  to  the  unique,  deterministic  jump  of  the  valence  electron  outside  the  closed  L- 
shell  in  Na  into  a  lower  orbit 

1  In  the  old  (Bohr)  quantum  theory,  an  atom  is  interpreted  as  radiating  like  an  oscillating  di]x>le,  and  from 
the  "Correspondence  Principle,"  fju  =  lu  fiu  =  1- 
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The  formula  (3.30)  is  discussed  in  Appendix  C,  which  presents  the  relation  between  the 
Einstein  A-  and  B-coefficients  for  spontaneous  and  stimulated  radiation  and  the  f-number. 

For  a  dipole-forbidden  transition  we  have 

l(ulril)|2  =  0  ,  (3.31) 

so  that  higher-order  terms  have  to  be  used  in  Eq.  (3.30b)  rather  than  (3.30a).  In  this  case, 
the  Aui  coefficient  is  evaluated  quantum  mechanically  for  all  lower  states  1  accessible  from  a 
given  excited  upper  state  u,  whose  lifetime  Tu  is  defined  as 

Tu-1/2iAui  .  (3.32) 

Now  we  can  define  an  effective  f-number  feff  for  the  forbidden  transitions  by  the  analog  to 
Eq.  (3.30b)  as 

feff.lu=  1.349  X  1021  Aui/v(u1)2  .  (3.33) 

For  a  dipole-forbidden  transition  of  the  type  considered  here,  the  standard  relation  (3.30a) 
between  the  dipole  matrix  element  and  the  f-number  does  not  hold,  although  formally 
Eq.  (3.30b)  can  always  be  used  to  define  an  effective  f-number  feff,  given  an 
A-coefficient.  Note  that  for  the  dipole-forbidden  transitions  of  Table  3-2B  the  effective 
transition  probabilities  Aui  or  feff  arc  much  smaller  than  for  allowed  transitions. 

If  there  are  ni  atoms  per  unit  volume,  each  of  absorption  cross  section  Oiu,  then  the 
absorption  coefficient  ka(v)  =  ni  aiu(v),  when  integrated  over  all  frequencies  contributing  to 
this  particular  transition,  gives  the  following  result.  The  line  strength  Siu,  which  in  the 
integral  of  the  absorption  coefficient  ka(v)  over  the  whole  of  a  particular  spectral  line  is 
defined  in  Appendix  C  by  the  relation 

Siu  = /ka(v)  dv  =  (7ce2/mc)  nabs  =  (tte2/mc)  ni  flu  >  (3.34a) 

where  v  is  frequency  (Hz).  In  other  words,  the  effective  number  of  absorbing  atoms  is 
Oabs  =  ni  fiu,  so  that  fiu  is  the  fraction  of  atoms  corresponding  to  absorption  in  the  1  =>  u 
transition.  We  can  of  course  write  Eq.  (3.34a)  in  the  equivalent  form 

i  oiu(v)  dv  =  nc  To  Jg(v)  dv  fiu  =  Ttc  ro  fiu  ,  (3.34b) 

where  ro  =  e2/mc2  =  2.82  x  10“^3  cm  =  "classical  radius  of  electron,"  and  g(v)  is  a 
normalized  spectral  shape  function  for  our  transition,  defined  so  that 

/  g(v)  dv  =  I  .  (3.34c) 

The  form  of  g(v)  under  different  conditions  will  be  discussed  below . 
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Table  3-2(a)  lists  transition  probabilities  for  some  dipole-allowed  transitions  of  H, 
N  and  O  atoms,  which  all  correspond  to  AL  =  ±1  (i.e.,  S  <->  P,  P  <->  D,  etc.,  but  not 
S  ->  S  or  S  ->  D  )  and  to  AS  =  0  (i.e.,  the  multiplicities  are  the  same  for  initial  and  final 
states).  Table  3-2(b)  lists  transition  probabilities  for  some  dipole-forbidden  transitions,  for 
which  either  AL  =  ±  2  (  S<->  D,  etc.)  or  AS  changes  (e.g.,  ^D-^P,  ^D-^S,  etc).  Note  that 
the  well-known  and  important  auroral  red  and  green  lines  of  atomic  oxygen  correspond  to 
dipole-forbidden  transitions.  The  point  is  that  the  allowed  transitions  of  N  and  O  are  either 
in  the  VUV  or  in  the  NIR,  while  these  forbidden  transitions  of  N  and  O  are  in  fact  quite 
important  because  of  the  large  column  densities  of  N  and  O  atoms  along  reasonable  lines  of 
sight,  although  their  transition  probabilities  are  much  smaller  than  those  of  optically 
allowed  transitions.'^^ 

Now  we  ask  for  the  character  of  the  shape  function  g(v),  which  is  typically  large 
near  the  line  center  and  falls  off  very  rapidly  to  either  side. 

There  are  three  principal  contributors  to  the  line  width: 

a.  Natural  line  width.  As  a  result  of  the  natural  lifetime  tn  of  the  upper  state  of 
the  transition  and  the  uncertainty  principle,  there  is  an  energy  width  AEn  = 
h/tN,  and  thus  a  natural  line  width 

ttN  =  AEN/h  =  l/27tTN  .  (3.35) 

b.  Collision  line  width.  The  upper  state  k  of  the  transition  is  deactivated  by 
collisions  with  other  particles  j,  with  number  density  nj,  collision  cross  section 
Oj,  and  velocity  uj 

ocoll  =  (l/27t)  Ij  nj  aj  uj  .  (3.36) 
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The  intensity  of  a  given  transition  is  characterized  either  in  terms  of  the  Einstein  A-coefficient  or  the 
oscillator  strength  or  f-number. 
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Table  3*2.  Transition  Probabilities  Aui.  Lifetimes  tu,  and  Absorption  Oscillator 
Strengths  f|u  for  Low-Lying  Atomic  States 

(a)  Allowed  Transitions  (Data  from  Radzig  and  Smirnov,  1985,  p.  225) 


Atom 

(ground  state) 

Transition 

(I-u) 

Wavelength 

X(nm) 

Aui  (10®secr‘') 

<lu 

gu/gi 

H(i2S) 

1S-2P 

121.57 

6.26 

0.416 

6/2 

3P 

102.57 

1.67 

6/2 

4P 

97.25 

0.681 

0.0290 

6/2 

2S-3P 

656.27 

0.225 

0.435 

6/2 

4P 

486.13 

0.0967 

0.103 

6/2 

2P-3S 

656.29 

0.0631 

0.0136 

2/6 

3D 

656.28 

0.647 

0.696 

10/6 

4S3 

486.14 

0.0258 

0.00305 

2/6 

4S 

486.13 

0.206 

0.122 

1 0/6  sic 

N(2p3-4S3/2) 

2‘*S3/2  •  3s^P 

120 

4.2 

0.27 

12/4 

2p^^P 

1.36 

0.080 

12/4 

2p3  2d  -  3s  2p 

3.4 

0.07 

6/10 

3S’  2d 

■9 

3.6 

0.08 

10/10 

0(2p^-3P2) 

2p4  3p  -  3s2S 

130.4 

6.0 

0.05 

3/9 

3d  3d 

102.7 

0.4 

0.01 

15/9 

3S  5S2  -3p  5p 

777.4 

1.0 

2.7 

15/5 

3S  3Si  -  3p  3p 

844.7 

0.28 

0.9 

9/3 

(b) 

Dipole-Forbidden  Transitions 

Atom 

Transition 

X(nm) 

Lite  Time  (sec) 

Aui  (sec-"*) 

Hi 

0(2p4) 

m 

557.7 

0.74 

1.35 

6.3  X  10-3 

630.0 

110 

9.09  X  10-3 

5.5x  10-''1 

■■ 

636.4 

N(2p3) 

2d -4s 

520.0 

9.4x104  (26  hrs) 

1.07  X  10-3 

4.3  X  10-'' 4 

2p-2D 

1040.0 

12 

8.33  X  10-2 

1.4  X  10-3 
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c.  Doppler  line  width.  Because  the  radiating  atom  moves  with  thermal  velocity 
(2kT/M)^/2^  there  is  ^  mean  width 

aD  =  (Vo/c)(2kT/M)l/2,  (3.37) 

where  v©  is  the  center  line  frequency,  T  the  temperature,  and  M  the  appropriate 
atomic  mass. 

Representative  values  for  cn,  ctcoll  and  ttD  as  a  function  of  altitude  in  the  earth's 
atmosphere  are  given  in  Fig.  3-4.  (These  calculations  were  actually  done  for  vibrational 
CO2  lines  in  the  4.3  pm  band,  rather  than  for  electronic  transitions.)  Note  that  on  is 
typically  very  small,  while  acoll  is  large  at  high  densities  and  predominates  in  the  earth's 
upper  atmosphere  at  altitudes  below  18  km;  at  higher  altitudes  (lower  densities)  Doppler 
broadening  determines  the  line  width.  Goody  (1964,  p.  99)  examines  line  shapes  in  more 
detail,  and  concludes  that  for  different  spectral  lines  the  altitude  above  which  Doppiei 
broadening  predominates  over  collision  broadening  ranges  from  7  km  to  37  km . 

These  different  mechanisms  lead  to  different  line  shape  function  g(v)  of  Eq.  (3.34). 
Natural  and  collision  broadening  lead  to  a  Lorentz  line  shape 

gL(v)  =  (ttL/ic)  /  [(v-Vo)2  -K  ttL^]  ,  (3.38a) 

where  ttL  =  ocn  +  a  Coll,  while  Doppler  broadening  leads  to  a  Gaussian  shape 

gD(v)  =  (l/aD7tl/2)  {exp -  [(v-Vo)2 /cd^]}  .  (3.38b) 

The  combination  of  Doppler  and  Lorentz  broadening  gives  rise  to  a  "Voigt  "line  shape— see 
e.g.,  Unsoeld,  1955,  p.  258f,  or  Arnold  et  al.,  1969,  esp.  p.  783. 

What  is  of  considerable  interest  is  how  the  effective  absorption  coefficient  changes 
with  optical  thickness,  the  so-called  Curve  of  Growth^^.  Let  us  define  a  mean  absorptance 
for  a  representative  spectral  line  in  a  band  by  the  relation 

Aav  =  (l/5)J[l-e-k(v.d]dv  ,  (3.39) 

where  5  =  mean  line  spacing  in  band  and  d  =  geometrical  thickness  of  radiating  medium. 


This  discussion  is  very  sketchy:  what  is  important  is  the  result  given  in  the-final  paragraph  of  this 
section.  A  more  detailed  discussion  is  given,  e.g.  in  Aller,  1963,  pp.  359-376. 

The  absorption  coefficient  A  must  not  be  confused  with  the  Einstein  A-coefficient 
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ALTITUDE  (km) 


Figure  3*4.  Natural,  Collision,  and  Doppler  Line  Widths  versus  Altitude 
(Calculations  Done  for  Lines  in  the  4.3  pm  Band  of  CO2) 
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We  now  ask  for  the  variation  of  Aav  with  the  integrated  optical  thickness^^ 

X  =  S  d  =  Jka(v)  d  dv  =  ni  flu  d  Jt  To  c  ,  (3.40) 

where  there  are  ni  atoms  per  unit  volume  in  the  lower  (absorbing)  state.  Figure  3-5  shows 
the  curve  of  growth,  which  is  the  general  character  of  the  absorption  coefficient  as  a 
function  of  integrated  optical  thickness  S  d/ao  and  of  the  ratio  2aL/otD  of  Lorentz 
broadening  to  Doppler  broadening.  Branch  QA  corresponds  to  the  optically  thin  case  for 
which 

Aav(opt.  thin)  ~  S  d  ~  ni  flu  d  .  (3.41a) 

Branch  QB  corresponds  to  the  Doppler  limit.  Because  of  the  Gaussian  line  shape  the 
absorption  saturates  rapidly  and  increasing  (S  d)  changes  Aav  very  little,  only 

Aav(Doppler)  “  ^log  (n,  f,^  d)  .  (3.41b) 

Branch  QC  corresponds  to  the  Lorentz  case.  Because  of  the  wider  wings  for  this  line 
shape,  Aav  increases  relatively  rapidly  with  S  d  or  ni  fiu  d 

Aav(Lorentz)  ~  (ni  flu  d)^^.  (3.41c) 

Note,  however,  that  the  limit  QD,  with  Aav  ~  (ni  flu  d)  for  large  optical  thickness,  is  not 
physically  realistic. 

The  discussion  should  be  referred  to  Section  3.3.1  [see  in  particular  Eqs.  (3.11) 
and  (3.12)]  for  the  two  limits  of  optically  thick  and  optically  thin  gas  parcels.  Typically, 
for  optically  thin  parcels  the  spectral  radiance  is  proportional  to  the  optical  thickness  or  to 
the  total  number  of  atoms  in  the  path,  while  for  optically  thick  parcels  the  spectral  radiance 
tends  to  a  black-body  radiance.  Note  that  the  total  absorption  and  radiance  due  to  a 
Doppler-broadened  line  increases  extremely  weakly  with  increasing  optical  depth  because  a 
Gaussian  profile  has  such  extremely  small  absorption  in  the  line  wings.  By  contrast,  a 
Lorentz  line  has  very  broad  wings  and  so  the  spectral  radiance  increases  far  more  rapidly 
with  increasing  optical  depth.  In  either  case,  the  blackbody  limit  is  evidently  attained  only 
asymptotically. 


'*5  The  optical  thickness  X  is  equal  to  kg'  x  of  Section  3.3.1. 


3.5  MOLECULAR  BAND  RADIATION 


In  contrast  to  atomic  spectral  lines,  which  correspond  to  the  transition  between  two 
electronic  states,  each  of  (fairly)  sharply  defined  energy  Eu  and  El  so  that  there  is  a  line  of 
ftequency 

v(ul)  =  (l/h)(Eu-Ei)  ,  (3.42) 

an  electronic  state  for  a  molecule  is  characterized  by  vibrational  quantum  number  v  and 
rotational  quantum  number  J  as  well  as  by  the  quantum  numbers  and  symmetry  factors  that 
characterize  the  electronic  state.  Because  the  nuclear  mass  is  so  much  larger  than  the 
electronic  mass,  the  motions  of  electrons  and  nuclei  can  be  decoupled— the  electrons  are 
treated  as  moving  in  the  field  of  temporarily  stationary  nuclei,  while  the  nuclei  move  in  the 
smeared-out  field  of  the  electrons.  This  is  called  the  Born-Oppenheimer  approximation. 
The  wave  function  is  written  as  the  product  of  electronic  and  vibrational/rotational  wave 
functions.  One  application  of  this  is  provided  by  Gilmore's  potential  energy  curves  for  the 
air  molecules  (see  Figs.  2-5  to  2-7).  The  energy  of  a  given  molecular  state  can  be  written 
as  the  sum  of  electronic,  vibrational,  and  rotational  contributions 

E(mol)  =  E(electronic)  +  coe  (v  -t- 1/2)  +  Be  J(J+1)  ,  (3.43) 

where  typically  E(electronic)  ~  1-5  eV  ,  coe  ~  2000  cm“^  (~  0.25  eV),  Be  ~  1  cm“^  to 
2  cm~l  (1.2  X  10^  -  2.4  x  10~^  eV  ).  The  vibrational  states  differ  by  ~  2000  cm"i  and 
have  rotational  structure  (AJ  =  0,1,2,...)  with  spacing  2  Be  ~  1-5  cm“l.  It  follows  from 
this  that  since  a  given  molecular  electronic  state  corresponds  to  numerous  vibrational  and 
rotational  levels,  a  transition  between  two  molecular  electronic  states  involves  very  many 
spectral  lines  characterized  by  the  vibrational  and  rotational  quantum  numbers  of  upper  and 
lower  electronic  states.  Typically,  the  population  of  highly  vibrationally  excited  states  will 
be  small,  but  the  population  of  highly  rotationally  excited  states  (with  rotational  quantum 
numbers  up  to  30)  will  be  large.  Thus,  we  recognize  that  a  molecular  electronic  transition 
is  not  characterized  by  a  single  spectral  line,  but  by  an  array  involving  thousands  of  spectral 
lines  of  varying  strength  and  width.  In  fact,  there  is  a  great  deal  of  structure  to  these 
arrays,  which  are  described  as  spectral  bands. 

Some  excellent  photographs  of  molecular  electronic  band  spectra  are  given  in 
Herzberg,  1950,  Ch.  n.  They  show  a  great  richness  of  phenomena  and  both  coarse  and 
fine  spectral "  ructure.  The  coarse,  banded  structtire  (~  2000  cm“l)  is  due  to  molecular 
vibrations,  while  the  fine  structure  due  to  molecular  rotational  lines  has  lines  spaced  1  cm~* 


to  5  cm~^  apart  In  general,  the  vibrational  structure  has  a  "sawtooth"  character,  falling  off 
rapidly  towards  one  side  (either  toward  higher  or  lower  frequency)  while  its  rate  of  decline 
toward  the  other  side  is  much  slower.  A  typical  electronic  band  system  covers  perhaps 
10,000-20,000  cm~l  (~  100  nm  in  the  visible),  so  that  it  has  ~  10,000  individual  lines,  of 
varying  spacing  and  strength.  From  the  discussion  of  Section  3.4  we  see  that  in  the  lower 
atmosphere  (below  10  km  to  40  km)  most  spectral  lines  are  collision  broadened  with  a 
typical  line  width  ~  0.01-0.1  cm“^  (see  Goody,  1964,  p.  99). 

In  the  IR  there  are  molecular  vibration-rotation  bands  due  to  molecules  such  as  NO, 
CO,  CO2,  H2O,  O3,  etc.  The  systems  typically  cover  1000  cm~l,  so  that  they  have 
perhaps  1000  or  less  important  lines  each,  but  because  of  the  differing  symmetry  of  these 
various  molecules,  and  more  particularly  because  of  the  applications  of  interest,  it  is  often 
desirable  or  necessary  to  treat  the  lines  in  real  detail. 

To  compute  the  transition  probability  for  a  molecular  transition,  one  again  uses  the 
Bom-Oppenheimer  separation  of  electronic  and  nuclear  motions,  in  which  the  molecular 
wave  function  is  described  as  the  product  of  an  electronic  and  a  vibrational  wave 
fiinction.^^  For  an  electronic  transition,  the  transition  probability  is  represented  in  terms  of 
the  transition  dipole  matrix  element  (ulrll)  where  r  =  electronic  coordinate  [see  Eq.  (3.30a)] 
as 

Transition  probability  -  fiu  qw’  .  (3.44) 

where  fiu  is  an  electronic  f-number  and  qwS  the  Franck-Condon  factor,  is  the  mean 
squared  overlap  of  vibrational  wave  functions  between  initial  and  final  states.  This  satisfies 
the  conditions 

qvv’  ~  2v'  qvv'  “  1  •  (3.45) 

If  the  mean  intemuclear  separation  of  a  molecule  in  its  initial  and  final  electronic 
states  is  very  similar,  the  vibrational  overlap  or  Franck-Condon  factor  qw'  is  largest  if 
Av  =  0,  with  transitions  with  Av  =  1,2,  having  successively  smaller  Franck-Condon 
factors.  However,  as  v  increases,  the  overlap  for  Av  =  0  decreases  and  so  does  the 
transition  probability  and  intensity  of  the  appropriate  spectral  lines.  All  this  structure  is 
implicit  in  the  energy  expression  of  Eq.  (3.43)  and  the  Bom-Oppenheimer  separation,  and 
it  gives  rise  to  the  structure  of  molecular  bands.  In  analogy  to  the  definition  of  an  atomic 


^  ®  From  the  correspondence  principle,  the  rotational  quantum  numbers  are  so  large  that  the  rotation  can  be 
treated  classically,  except  for  H2,  HD,  and  D2  at  low  temperatures. 
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line  strength  Siu  of  Eq.  (3.34a),  the  transition  probability  for  a  band  system  is  normally 
described  in  terms  of  the  integrated  band  intensity  Siu 

Siu  =  /k(v)dv(cm-2atm-l)  ,  (3.46) 

or  equivalently,  a  band  f-number,  which  is  related  to  the  band  strength  by  the  relation 

Siu  =  2.3795  X  107  flu  •  (3.47) 

Appendix  C  presents  the  numerical  factors  in  terms  of  fundamental  constants  and  ties  the 
various  parameters  together  in  a  schematic  way. 

As  in  the  case  of  an  atomic  transition,  there  are  selection  rules  for  molecular 
transitions,  but  these  are  naturally  more  complex  because  there  are  so  many  more  quantum 
numbers  and  other  symmetries  involved. 

For  electronic  transitions  of  diatomic  molecules,  we  have  the  selection  mles'‘7; 


AA  =  0,  ±  1  (E  ->  E;  1  — >  11;  11  ->  D,  but  not  E  ->  A)  (3.48) 

AS  =  0  (nA  ->  n'A’  only  if  n'  =  n)  (3.49) 

+  <=>  -  ,  not  +=>  +  ;-=>-  (3.50) 

g  <=>  u  ,  not  g=>  g ,  u  =>  u  .  (3.51) 


For  vibrational  transitions  within  a  given  electronic  state,  things  are  quite  different. 
A  homogeneous  diatomic  molecule  like  N2  or  O2  has  no  dipole  moment,  so  that  there  are 
no  transitions  between  different  vibrational  levels.  For  a  heteronuclear  diatomic  molecule 
like  NO,  the  selection  rule  Av  =  ±  1  is  only  weakly  obeyed  because  molecules  do  not 
behave  completely  like  harmonic  oscillators  for  which  this  selection  rule  holds  rigorously. 

For  rotational  transitions  in  homonuclear  diatomic  molecules,  the  selection  rule 
AJ  =  ±  1  is  strictly  obeyed.'** 

Table  3-3(a)  lists  characteristics  for  some  important  electronic  transitions  in 
diatomic  air  molecules  which  are  mainly  in  the  UV  and  visible,  and  Table  3-3(b)  lists 
similar  characteristics  for  vibrational  transitions,  which  are  always  in  the  IR.  Table  3-3(a) 
presents  a  minimal  characterization  of  six  electronic  band  system,  used  for  an  approximate 


'*7  Sec  Appendix  A  for  a  discussion  of  the  various  quantum  numbers  and  other  symmetry,  parameters. 

Pure  rotational  transition*!  are.  nortnaHy  in  the  millimeter  wave  spectral  region  which  is  not  normally 
of  interest  here,  but  rotational  structure  occurs  in  vibrational  and  electronic  band  systems. 


Table  3-3.  Some  Characteristics  of  Important  Molecular  Band  Systems  of  Air^^ 

(a)  Electronic  Transitions 


Band  System 

Wavelength  ( ttm) 

Sband(cm-2atm-‘') 

flu 

N2  (First  Positive, B3n->A3l ) 

0.52-1.1 

6.7  X  10^ 

0.0028 

N2  (Second  Positive,  C^n->B^n) 

0.30-0.47 

1.0  X  10® 

0.043 

N2'^  (First  Negative,  b2i->x2i) 

0.34-0.52 

8.3  X  10® 

0.035 

O2  (Schumann-Runge, b2z->x2z)  - 

0.19-0.43 

(6x10®) 

0.17-0.34X. 

NO  (Beta  ,B2n->x2n) 

0.21-0.52 

1.3  X  10® 

0.0053 

NO  (Gamma>\2x.>x2n) 

0.19-0.28 

5.9  X  10^ 

0.0025 

Note: 

a.  the  wavelength  range  is  representative  only 

b.  Sband  (which  is  sometimes  called  Siu)has  some  temperature  dependence. 

c.  flu  is  sometimes  described  as  absorption  f-number,  fabs>  to  distinguish  it  from  the 
emission  f-number  fyi  [see  Appendix  C,  Eq.  (C.f  0)]. 

d.  Note  that  the  f-number  for  the  O2  Schumann-Runge  band  has  a  strong  wavelength 
dependence  (X-values  in  p.m): 

f|u  »  0.10  at  X.-  0.2  |im 

f|u  »  0.034  at  X. »  0.4  |im. 


These  data  are  compiled  from  various  sources  and  should  be  considered  as  representative  only.  For 
more  up-to-date  values  of  electronic  transitions  (Part  (a))  see  Kuz’menko  et  al,  1979;  and  for  vibrational 
transitions  (Part  (b)),  see  Kennealy  &  Del  Greco,  1983. 
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Table  3*3.  (continued) 


(b)  Vibrational  Transitions 


Molecule 

v->v' 

Wavelength  (pm) 

Sband  (cm“2  atm""') 

flu 

CO 

0->  1 

4.67 

260 

1.09  X  10-5 

0-5.2 

2.34 

1.79 

7.5x10-® 

NO 

5.33 

110 

4.7  X  10-6 

2.69 

1.3 

5.4  X  10-8 

OH 

2.8 

100 

4.0  X  10-6 

BBI 

1.4 

4.0 

1.7  X  10-7 

CO2 

ooo->  010 

15 

187 

7.9  X  10-6 

ooo->  001 

4.3 

2700 

1.14  X  10-4 

000->  101 

2.69 

39 

1.6  X  10-6 

000->  201 

2.0 

1.0 

4.0  X  10-8 

H2O 

000->  010 

6.27 

300 

1.3  X  10-6 

000->  001 

2.66 

100 

4.5  X  10-6 

000->  01 1 

1.88 

30 

1.3  X  10-6 

000->  101 

1.38 

20 

8.0  X  10-7 

000->  1 1 1 

1.14 

0.6 

2.5  X  10-8 

O3 

000->  001 

9.6 

260 

1.1  X  10-6 
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description  of  the  radiation  from  high-temperature  air.  These  (electronic)  f-numbers  are 
quite  small  (10~2-1(>-1 )  compared  to  the  f-numbers  for  the  Na-D  or  H-Ly-a  lines,  but  they 
are  in  the  general  range  of  the  allowed  atomic  transitions  shown  in  Table  3-2(a).5®  Note 
that  while  five  of  the  six  band  systems  in  Table  3-3(a)  are  relatively  well-behaved,  the  O2 
Schumann-Runge  system  is  anomalous,  covering  a  very  wide  spectral  range  and  having 
both  a  wavelength-dependent  f-number  and  a  very  non-diagonal  Franck-Condon  array. 
The  reason  for  this  is  that  the  intemuclear  separation  of  the  O2  B  and  X  states  is  quite 
different  (1.60  and  1.21  A  respectively)-see  Herzberg,  1950,  p.  560,  and  also  Fig.  2-7. 

Because  of  their  symmetry,  homonuclear  diatomic  molecules  like  N2  and  O2  do  not 
have  IR  active  vibration-rotation  transitions,^^  but  NO  and  the  triatomic  species  CO2,  H2O, 
O3,  N2O,  and  NO2  are  important  radiators  in  the  infrared  which  are  indicated  in 
Table  3-3(b).  Note  that  the  transition  probabilities  (f-numbers)  for  the  vibration-rotation 
bands  are  significantly  smaller  than  those  for  electronic  bands. 

How  the  band  structure  of  a  given  molecular  transition  in  a  specified  spectral  range 
is  treated  depends  on  the  particular  application.  Thus,  if  one  asks  for  the  emission  from 
high-temperature  gases  at  low  spectral  resolution,  very  little  detail  is  necessary.  In 
particular,  for  an  optically  thin  situation  when  self-absorption  and  scattering  are 
unimponant  so  that  a  Planck-type  mean  (see  Eq.  (3.27)]  is  all  that  one  needs,  the  total 
integrated  intensity  suffices.  Thus,  the  description  in  terms  of  electronic  system  f-numbers 
with  vibrational  structure  analyzed  in  terms  of  a  Franck-Condon  factor  is  satisfactory  for 
determining  the  total  radiation  from  air  or  combustion  products  at  high  temperatures. 

However,  if  one  asks  for  the  atmospheric  transmittance  through  1  km  to  10  km  path 
lengths  (or  through  the  whole  atmosphere  in  vertical  viewing),  especially  at  high  spectral 
resolution  in  the  IR,  where  there  are  numerous  absorption  regions  alternating  with 
"windows,"  much  more  detail  is  needed.  One  example  of  this  is  when  one  views  not  a 
gray  body  through  a  long  atmospheric  path,  but  rather  a  rocket  plume  or  other  array  of  hot 
CO2  and  H2O  gas  viewed  through  the  atmosphere  containing  cold  H2O  and  CO2  gas,  or 
for  laser  transmission  at  very  high  spectral  resolution. 


The  reason  why  these  f-numbers  for  given  transitions  are  small  is  that  the  upper  states  radiate  into  a 
large  number  of  lower-lying  states  over  a  wide  spectral  range:  it  is  unusual  for  upper  states  such  as 
Na  (3^P)  to  radiate  into  just  one  lower  state  (here  3^S). 

^  ^  But  they  do  have  Raman  spectra,  which,  however,  are  not  discussed  here. 
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During  the  1950s  and  1960s  a  number  of  Band  Models  were  devised  and  tested 
against  experiment.  These  models  use  lines  of  equal  or  random  spacing,  of  uniform  or 
varying  strength,  sometimes  with  a  configuration  like  the  "Q-Branch"  of  a  band  [see,  e.g., 
Herzberg  (1950),  p.  48,  etc.].  Two  limits  are  the  weak-  and  strong-line  approximations. 
In  the  weak-line  limit  the  line  is  optically  thin,  even  at  the  line  center,  so  that  the  emission 
or  absorption  is  proportional  to  the  number  of  molecules  in  the  path  [Aav  ~  (Nfd)"See 
Eq.  (3.41a)  or  branch  AQ  of  Fig.  3-5].  By  contrast,  in  the  strong-line  limit  the  line  is 
optically  thick  at  the  center,  so  that  increasing  the  number  of  molecules  in  the  path  simply 
increases  the  emission  or  absorption  in  the  line  wings,  and  of  course  this  gives  a  different 
variation  of  absorption  Aav  with  Nfd,  depending  on  the  line  shape  (Doppler,  Lorentz,  or 
Voigt,  see  Fig.  3-5);  the  strong-line  approximation  traditionally  corresponds  to  the  Lorentz- 
line  limit  Aav  ~  (Nfd)l/2  [^f.  Eq.  (3.41c),  or  branch  QC  in  Fig.  3.5].  For  a  review  of  band 
models,  see,  e.g..  Goody,  1964,  Ludwig  et  al.,  1973. 

With  the  great  increase  in  computer  power  over  the  last  20  years,  it  has  become 
feasible  to  describe  vibration-rotation  bands  in  the  IR  by  a  line-by-line  analysis,  rather  than 
by  using  band  models.  Note,  however 

(a)  For  electronic  transitions,  it  is  still  frequently  convenient  to  use  band  models 
rather  than  line-by-line  descriptions,  e.g.,  in  the  AFGL  FASCODE 
transmission  code,  because  of  the  very  large  number  of  spectral  lines. 

(b)  For  many  applications,  it  is  convenient  to  describe  vibration-rotation  bands  by 
band  models  rather  than  by  line-by-line  analyses,  partly  because  it  is  quicker 
and  partly  because  the  detailed  input  data  needed  for  line-by-line  treatments  is 
not  available.  Thus,  in  the  1970s  an  Atlas  of  Spectral  Lines  was  developed  at 
AFGL.  By  now  this  Atlas  (Rothman  et  al.,  1987)  contains  in  excess  of 
350,000  individual  lines,  and  in  the  IR  line-by-line  calculations  can  now  be 
carried  through  for  computations  which  could  not  have  been  considered 
20  years  ago,  by  using  the  FASCODE  computer  program.  However,  this 
HITRAN  listing  is  quantitatively  correct  only  up  to  -  600  K;  at  higher 
temperatures  there  are  many  more  spectral  lines  of  importance.  Thus,  while 
FASCODE  may  give  reasonable-appearing  results  for  high-temperature  air  and 
other  species,  yet  some  important  collections  of  lines  will  be  missing.  The 
HITEMP  compilation  consisting  of  CO2  and  H2O  lines  in  selected  absorption 
bands  for  temperatures  up  to  1000  K  to  1500  K  is  being  prepared  at  AFGL, 
but  has  not  yet  been  published. 
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Let  us  now  present  some  figures  for  the  absorption  cross  sections  of  various 
molecules.  Figure  3-6  shows  the  O2  absorption  coefficient  in  the  far  UV  spectrum.  O2  can 
dissociate  into  two  ground  state  (^P)  oxygen  atoms: 

02(X3S)  =>  20(3P)  ,  (3.52a) 

for  X  <  242  nm  (5.1  eV  =  O2  dissociatioii  energy)  and  the  weak  Herzberg  continuum  is 
found  there.  The  Schumann-Runge  (SR)  bands  fall  generally  from  175  nm  to  195  nm: 
below  175  nm  (7.1  eV)  there  is  the  strong  Schumann-Runge  continuum: 

02(X3Z)  =>  0(3P)  +  0(lD)  (7.1  eV)  ,  (3.52b) 

and  below  ~  130  nm  the  02(X3X)  molecule  can  dissociate  into  0(3p)  -1-  0(^S).  Figure  3-6 
shows  the  short  wavelength  portion  of  the  SR  band  system  (for  X  >175  nm)  and  the  SR 
continuum  (140  nm  to  175  nm).  For  X  <  125  nm  there  is  evidently  at  least  one  more  band 
system,  but  this  is  not  important  for  our  purposes.  At  longer  wavelengths  the  cross 
sections  become  smaller  with  a  threshold  (  for  the  Schumann-Runge  B  ->  X  transition)  at 
242  nm. 

Figure  3-7  shows  die  ozone  spectral  absorption  in  the  visible  and  the  near  UV.  The 
strongest  system  is  the  Hartley  system  which  includes  the  photodissociation  of  the  non¬ 
linear  ozone  molecule 

03(XlAi)=>02(alA)  +  0(lD)  (3.99  eV,  X  <  31 1  nm)  .  (3.53) 

which  is  very  important  for  the  absorption  of  solar  UV  in  the  earth’s  atmosphere. 
Figure  3-7(a)  shows  the  cross  section  of  O3  in  the  near-UV;  the  strong  Hartley  system  lies 
below  X  =  300  nm,  while  the  longer-wave  tail  (i.e.,  X  >  300  nm)  is  known  as  the  (weak) 
Huggins  system  (C-X).  Ozone  also  absorbs  in  the  orange/red  portion  of  the  visible:  the 
weak  Chappuis  (B-X)  absorption  cross  section  is  shown  in  Fig.  3-7(b). 

By  comparing  the  maximum  absorption  cross  section  in  Figs.  3-6  and  3-7,  we  see 
that  a  strong  absorption  corresponds  to  a  cross  section  of  10~^*  to  lO"!’^  cm^,  while  a 
weak  absorption  is  down  by  perhaps  a  factor  10  to  100. 

3.6  ATOMIC  RADIATION  IN  GENERAL;  FREE-FREE,  FREE-BOUND, 
BOUND-BOUND 

We  have  given  a  rather  careful  and  detailed  discussion  of  the  effective  spectral 
absorption  arising  from  a  single  atomic  spectral  line,  and  from  the  electronic  band  system 
of  a  diatofnic  molecule,  which  is  an  array  of  separate  lin>.  Individual  atomic  lines  are  not 
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Figure  3-6.  O2  Absorption  Cross  Section  in  the  UV. 
(Source:  Watanabe,  1958) 
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Figure  3-7.  O3  Absorption  Cross  Section; 

(a)  UV  (conventionally  one  speaks  of  the  Hartley  system  below  300  nm 
and  the  Huggins  system  at  longer  wavelengths); 

(b)  Chappuis  system  In  the  Visible. 

(Source:  VIgroux,  1953"See  also  Goody,  1964) 
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particularly  important  from  the  standpoint  of  energy  transfer  involving  blackbody  radiation 
because  a  line  is  so  narrow  that  the  total  energy  carried  within  its  width  is  a  small  fraction 
of  the  total  Planck  spectrum,  and,  further,  it  takes  very  few  atoms  to  saturate  a  spectral 
line.^^  However,  at  temperatures  above  6000  K,  the  N  and  O  (and  other)  atoms  have  a 
sufficiently  large  number  of  lines  and  continua  that  their  radiation  is  important. 

From  the  standpoint  of  the  overall  thermal  radiation  from  high-temperature  air,  one 
is  most  interested  in  absorbers  covering  a  wide  spectral  interval,  either  arrays  of  individual 
spectral  lines  (such  as  molecular  band  spectra,  which  have  been  discussed  in  Section  3.5) 
or  in  continua,  which  are  discussed  here.  The  cross  sections  involved  here  are  typically 
very  small,  but  the  large  available  column  densities  of  potential  absorbers  more  than  make 
up  for  the  narrow  spectral  interval  covered  by  atomic  resonance  lines. 

A  typical  atomic  energy  level  diagram  is  shown  in  Fig.  3-8.  Note  that  there  are  a 
number  of  discrete  levels  0,1, 2,3,...  which  cluster  towards  an  ionization  limit  J  and  which 
make  up  the  discrete  spectrum  D:  transitions  "a"  and  "b"  belong  to  the  discrete  spectrum, 
and  may  be  called  bound-bound  transitions.  Above  the  ionization  limit  J  all  energy  values 
are  allowed,  giving  the  continuous  spectrum  Transitions  involving  at  least  one 

state  in  the  continuum  are  not  characterized  by  sharp  spectral  lines,  because  all  energy 
values  in  the  continuum  are  available.  We  distinguish  between  transitions  involving  one 
bound  state  and  one  free  state,  such  as  "c",  which  is  a  free-bound  transition,  and  transitions 
whose  initial  and  final  states  are  both  in  the  continuum,  such  as  the  free-free  transition  "d". 

What  is  the  physical  meaning  of  these  transitions?  The  free-bound  transition  "c" 
corresponds  to  a  radiative  ionization/deionization  process: 

A(l)-i-hv=>  A+-^e-  ,  (3.54) 


For  instance,  consider  the  Na-D  line,  actually  a  doublet  at  589.0-589.6  nm.  At  STP,  the  collision  line 
width  is  given  by  Eq.  (3.36)  as  ocoll  ~  10^  Hz,  and  thus  from  Eq.  (3.34b)  we  have  Oabs  ~  10“' '  cm^ 
within  a  line  of  0.001  nm  of  both  589.0  and  589.6  nm.  Now,  if  the  number  density  of  Na  is  just 
times  ambient,  then  the  line  is  optically  thick  already  in  a  path  length  of  4  cm.  Molecular  and 

continuum  cross  sections  are  in  the  range  10“^®  to  10“'®  cm^,  but  this  is  clearly  quite  sufficient  for 
the  line  center  to  be  optically  thick  for  a  path  going  through  a  significant  fraction  of  the  atmosphere. 

Notice  the  narrow  region  K  Just  below  the  ionization  limit  J.  Here  there  are  closely-spaced  discrete 
(i.e.,  bound)  levels,  but  at  sufficiently  large  density  or  high  temperature  where  the  levels  are  broadened 
sufficiently  so  that  the  level  width  is  greater  than  the  level  spacing  gives  rise  to  a  quasi-continuous 
spectrum  or  an  effective  reduction  in  the  ionization  potential. 

^  There  are  also  auto-ionizing  discrete  states  which  have  structure  above  the  ionization  limit 
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while  a  free-free  transition  corresponds  to  the  interaction  of  a  free  electron  with  an  A+  ion, 
in  which  the  electron  is  decelerated  in  the  field  of  the  ion,  radiating  its  excess  energy  as  a 
light  quantum: 

A+ +  e- +  KE  =>  A+ -I- e- -I- hv -I- KE'  .  (3.55) 

This  Bremsstrahlmg  is  of  considerable  importance  in  the  context  of  air  opacity,  especially 
since  the  effect  is  observed  also  for  an  electron  in  the  field  of  a  neutral  atom  with  which  it 
interacts  by  polarization  or  other  shorter-range  interactions. 

In  practice,  one  treats  all  these  continuum  terms  together— the  free-free,  free-bound 
as  well  as  the  closely  spaced  lines  corresponding  to  the  series  limit  region  K,  and  for  the 
case  of  a  hydrogen-like  atom  (which  has  one  electron  outside  a  nucleus  of  charge  +Ze)  the 
expression  for  the  total  radiance  L(A,,T)  due  to  all  these  three  terms  has  a  simple  form: 

L(X,T)  =  C*  Z2  ne  ni/(kT)l/2  ;  C*  =  6.36  x  10^7  cgs  ,  (3.56) 

for  L  (X,T)  of  Eq.  (3.56)  in  appropriate  cgs  (not  MKS)  units.^s  For  a  singly  ionized  atom, 
the  effective  charge  Z  =  1;  and  for  the  neutral  free-free,  etc.,  continuum,  such  as  the 
reaction 

A  e“  +  KE  =>  A  +  e”  -t-  hv  +  KE'  ,  (3.57) 

one  can  still  use  the  expression  (3.56),  but  with  an  appropriate  numerical  value 

Zeff2  =  Zeff2(X,'i:)  .  (3.58) 

Thus  for  T  ~  10,000  K  in  air,  Zeff^  -  1  in  the  visible  and  ~  0.1  in  the  near  infrared.^^  For 
the  case  of  electron-neutral  Bremsstrahlmg  one  puts  n^  ntot  in  Eq.  (3.56)  instead  of  n©  nj 
(where  ni  =  ion  number  density,  ntot  =  number  density),  and  thus  for  a  sufficiently 
small  degree  of  ionization,  this  neutral  free-free,  etc.,  radiation  is  actually  more  important 
than  that  from  the  ionized  continuum.57 

3.7  ABSORPTION  COEFFICIENT  OF  IIIGH-TEMPERATURE  AIR 

The  absorption  coefficient  of  high-temperature  air  has  often  been  treated  as  due  to 
the  canonical  six  electronic  band  systems  listed  in  Table  3-3(a),  plus  the  NO  vibration- 


^  ^  Here  the  number  densities  n^  and  nj  are  listed  per  cm^. 

See,  for  example,  Morris  et  al.,  1969,  and  Zeldovich  and  Raizer,  1966,  Vol.  I,  p.  255  ff. 

57  This  explains  why  the  H  to  H“  continuum  predominates  over  H''"  to  H  in  the  solar  photospheric 
emission. 
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rotation  bands  and  a  number  of  continua.  Table  3-4  lists  the  continuous  contributors  that 
are  important  at  relative  density  p/po  =  10"^  in  the  computations  of  Landshoff  and  Magee, 
1969  (who  also  include  the  N2  Birge-Hopfield  b'lZu'^  -  XlZg+  electronic  band  system  in 
their  calculations). 


Table  3-4.  Continuous  Contributors  for  the  Absorption  Coefficient  of 
High-Temperature  Air  (Source:  Landshoff  and  Magee,  1969) 


System 

Spectral  Region  of  Importance 
X(nm)  Energy  (eV) 

Temperature  Region 
(K) 

O2  Schumann-Runge  Continuum 

130-175  7.1-9.5 

1000-10000 

NO2 

240-1130  -1. 1-5.1 

1000-5000 

O'  pnotodetachment 

120-830  -1.5-10.7 

>  3000 

N  pfotoionization 

Entire  range  of  calculation 

>  7000 

0  piiotoionization 

Entire  range  of  calculation 

>  7000 

Fret-free  in  presence  of  ions 

>  1600  <  0.8 

Wider  range  at  higher 
temperatures 

>  6000 

The  computation  is  quite  straightforward,  if  laborious.  First  one  determines  the 
equil-brium  chemical  composition  of  the  material,  then  identifies  the  appropriate  radiating 
syste  IS  and  computes  the  spectral  absorption  coefficient  at  moderately  low  resolution, 
usin’  an  appropriate  band  model.  In  Fig.3-9  we  show  the  spectral  absorption  coefficient 
for  V  jious  density  ratios  D  =  p/po  at  temperatures  T  =  2000,  4000,  6000,  9000,  and 
1200c  K  from  the  results  listed  in  Landshoff  and  Magee,  1969. 

The  following  points  should  be  noted: 

1 .  In  an  optically  thin  medium,  the  absorption  coefficient  kx.  is  proportional  to  p 
or  p,  as  a  measure  of  the  number  of  absorbing  or  emitting  molecules. 

2.  The  temperature  dependence  is  largely  determined  by  the  species  that  radiate. 
At  intermediate  temperatures  there  are  more  metastable  molecules  (such  as  NO, 
CN,  OH,  etc.)  which  tend  to  be  effective  radiators  than  either  at  low 
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Waveiangth  (um) 

(e)  T  r  12,000  K 
Figure  3-9.  (continued) 

temperatures  (below  -  1000  K),  where  there  are  few  reactive  molecules,  or  at 
high  temperatures  where  there  are  few  molecules  left.  (See,  e.g..  Fig.  2-3  in 
Chapter  2). 

At  the  highest  temperatures  all  the  molecules  dissociate,  giving  rise  to  a 
spectrally  rather  uniform  absorption  coefficient  due  to  free-free  and  free-bound 
continua.  [See,  e.g..  Fig.  3-9(e):  at  12,0(X)  K  there  are  very  few  molecules 
left.] 

In  addition  to  various  continua,  a  large  number  of  molecular  band  systems 
contribute  over  the  temperature  range  considered  here.  Thus,  the  later 
computations  of  Johnston  et  al.,  1977,  use  21  electronic  band  systems  in 
addition  to  the  NO  vibration  and  a  number  of  continua,  including  those  of 
Table  3-4  and  some  others.  Note,  however,  that  older  air  calculations  using 
just  6  band  systems  agree  fairly  well  (factor  of  two)  with  the  observations, 
because  the  f-numbers  are  obtained  from  measurements  of  the  radiation  from 
high-temperature  air. 
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3.8  RADIATION  FROM  HIGH-TEMPERATURE  AIR 


Figure  3-10  gives  the  energy  radiated  from  hot  air  as  a  function  of  temperature  and 
density.  For  ICBM  or  spacecraft  reentry,  peak  heating  typically  occurs  in  the  altitude  range 
15  km  to  50  km,  or  p/po  ~  0.16  to  0.00084,  so  that  the  p/p©  values  10“3  to  10"^  are  most 
relevant  to  flight  conditions.  Kivel,  1961,  points  out  that  radiative  heating  will  be 
important  as  compared  to  aerodynamic  heating  for  high  velocities  (>  9  km/sec)  and  low 
altitudes  (<  45  km).  For  the  case  of  vehicles  with  ablative  heat  shields,  radiation  may 
become  important  at  lower  velocities  and  possibly  at  higher  altitudes  because  of  the 
production  of  small  particles  (e.g.,  soot)  which  are  much  more  effective  radiators  than  air 
molecules.58  In  Fig.  3-10  we  also  show  the  black-body  radiation  per  unit  area  (W/cnfi) 
which  is,  of  course,  only  comparable  with  radiative  gas  heating  in  W/cm^  if  we  refer  to  a  1- 
cm  thick  layer  of  air.  Since  the  shock  layer  (or  boundary  layer)  is  roughly  1  cm  thick,  the 
comparison  is  not  inappropriate,  but  can  in  any  case  be  used  for  scaling  purposes.  Since 
stagnation  temperatures  for  satellite  or  ICBM  reentry  lie  in  the  range  6000  K  to  8000  K,  we 
see  that,  indeed,  radiative  heating  due  to  air  is  quite  small  compared  to  blackbody  radiation 
from  the  surface. 


Ablative  heat  shields  such  as  are  used  in  many  reentry  vehicles  can  have  quite  different  radiative 
properties  from  pure  air.  Thus,  at  1  atmosphere  pressure  and  5000  K,  the  Planck  mean  absorption 
coefficient  of  a  gaseous  50  percent  carbon/50  percent  air  equilibrium  mixture  is  2500  times  as  large  as 
that  of  pure  air,  while  if  the  carbon  is  present  as  small  soot  particles,  the  Planck  mean  is  larger  still  by 
another  factor  of  15  (see  Main  and  Bauer,  1966,  1967).  Therefore,  a  1-cm  thick  soot  layer  would 
radiate  like  a  blackbody. 
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4.  SOME  HEATING  MECHANISMS 


4.1  INTRODUCTION 

Chapter  2  examines  the  composition  of  high-temperature  air  in  Local 
Thermodynamic  Equilibrium  (LTE),  and  Chapter  3  considers  the  radiation  from  high- 
temperature  air,  in  effect  in  Radiative  Steady  State  (RSS).  These  are  limits  which  may  not 
be  attained  in  practice,  but  which  are  useful  for  specific  applications.  Here  we  present 
several  different  applications.  For  each  one  we  ask  how  energy  is  transferred  into  the 
system  and  hence  how  adequately  or  rapidly  is  LTE  (or  RSS)  approached. 

The  following  topics  are  addressed  here: 

•  Reentry  heating,  which  is  due  to  aerodynamic  (shock  and  boundary  layer) 
heating,  is  treated  in  Section  4.2. 

•  Nuclear  fireball:  this  entails  heating  by  keV  X-rays,  and  is  discussed  in 
Section  4.3. 

•  Upper  atmosphere-photodissociation  of  O2  and  O3  by  solar  UV  is  discussed  in 
Section  4.4.  Unlike  the  nuclear  and  reentry  cases,  this  is  clearly  a  non-LTE 
situation 

Light  propagation  effects  are  not  considered  here,  but  they  can  involve  processes  which  are 
very  far  from  LTE,  in  particular  for  the  case  of  high-energy  lasers. 

All  these  examples  demonstrate  how  the  energy  of  excitation  moves  towards  LTE. 
The  details  of  the  equilibration  process  are  highly  specific  for  each  particular  example, 
proceeding  through  some  particular  inelastic  collisions  and  rate  processes.  It  is  generally 
appropriate  to  introduce  a  characteristic  time  for  energy  transfer 

"^har  ~  l/(^)  (4.1) 

*  where  n  (cm-3)  is  the  number  density  of  colliding  particles  (~  2.5  x  1 0^9  cm“3  at 

1  atmosphere  pressure— see  Table  2-4  on  p.  2-22),  while  the  rate  constant  k  (cm^/sec) 
refers  to  a  specific  rate-limiting  step  and  is  the  product  of  the  reaction  cross  section  a  (cm2) 
and  of  the  molecular  collisional  velocity  v  (cm/sec).  Typically  a  ~  10~20 .  io-l5  cin2  and 
I  V  ~  105  cm/sec,  so  that  k  =  ctv  ~  lO'^^ .  jq-IO  cm^/sec. 
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If  the  critical  species  is  present  to  1  part  per  million,  so  that  n  ~  10^^  cm~3,  then  the 
charactenstic  time  tchar  is  in  the  range  1  msec- 100  sec,  but  the  actual  range  if  often  much 
larger.  In  general,  the  system  attains  LTE  in  times  much  greater  than  Xchar- 

4.2  REENTRY  HEATING 

An  ICBM  or  low-earth-orbiting  satellite  travels  in  space  at  a  speed  on  the  order  of 
7  km/sec.  As  the  vehicle  enters  the  earth’s  atmosphere,  it  is  slowed  down  by  collision  with 
air  molecules  whose  mean  thermal  speed  is  ~  10^  cm/sec  =  1  km/sec.  This  leads  both  to  a 
hypersonic  shock  wave^9  that  stands  off  in  front  of  the  vehicle  and  also  a  boundary  layer^o 
on  the  vehicle.  Figure  4-1  sketches  how  the  initial  kinetic  energy  of  motion  of  the  vehicle 
is  converted  into  heat.  Figure  4- 1(a)  shows  schematically  how  the  atmospheric  density  p 
and  the  vehicle  velocity  u,  deceleration^ '  and  heating  rate^^  change  with  altitude,  and 
Fig.  4- 1(b)  presents  the  energy  conversion  from  kinetic  energy  into  heat  in  semi- 
quantitative  form.  It  shows  that  fraction  of  the  kinetic  energy  of  the  vehicle  which  reaches 
the  surface  as  heat.  As  the  vehicle  comes  down  below  100  km  where  the  gas-kinetic  mean 
free  path  is  16  cm,  a  bow  shock  gradually  builds  up,^^  and  through  this  shock  (and  the 
boundary  layer  which  develops  on  the  surface  of  the  vehicle)  the  kinetic  energy  of  the 
vehicle  is  gradually  convened  into  heat,  some  of  which  reaches  the  surface.  Note  that  as 
the  vehicle  slows  down,  its  energy  conversion  fraction  decreases.  At  altitudes  in  the 
30-km-altitude  range,  the  flow  along  the  body  changes  from  laminar  to  turbulent,  with  a 
significantly  greater  heat  transfer  rate.  For  reentry  velocities  of  approximately  7  km/sec,^ 
radiative  heating  is  low  because  at  temperatiu’es  below  -  12,000  K  the  radiation  from  hot 


When  a  vehicle  travels  through  the  atmosphere  faster  than  the  local  sound  speed  (which  is 
~  0.3  km/sec  at  100  km  altitude),  a  shock  wave  develops  because  the  atmosphere  cannot  accommodate 
to  the  motion  in  a  less  drastic  fashion.  As  the  vehicle  speed  increases,  the  shock  wave,  which  initially 
is  almost  perpendicular  to  the  direction  of  motion  of  the  vehicle,  turns  so  as  to  make  successively 
smaller  angles  to  the  velocity  vector.  At  hypersonic  conditions  the  shock  '  wraps  around"  the  front  of 
the  body. 

The  boundary  layer  is  a  region  of  a  supersonic  flow  close  to  a  solid  surface  in  which  most  of  the 
interaction  between  the  surface  and  the  flow  occurs.  Thus  the  flow  here  is  subsonic  and  non-ideal  gas 
effects  such  as  viscous  dissipation  and  ablation  occur  in  the  boundary  layer. 

^  ^  The  deceleration  is  proportional  to  pu^. 

^2  The  heating  rate  is  proportional  to  pu^. 

At  higher  altitudes  where  the  vehicle  size  is  smaller  than  the  mean  free  path,  the  presence  of  the  vehicle 
does  not  affect  the  air  flow  and  so  no  shock  develops. 

This  corresponds  to  reentry  from  low  earth  orbit  or  to  ICBM  entry. 
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Figure  4-1,  Energy  Transfer  on  Reentry. 

(a)  Schematic  of  Slowing-Down  and  Heating-Up  as  Function  of  Altitude, 
(b)  Energy  Conversion  for  Reentry  at  Different  Velocities. 


air  is  small— see  Fig.  3-10  on  p.  3-41.  However,  at  higher  entry  velocities— such  as  11 
km/sec,  as  is  illustrated  in  Fig.  4-l(b)-radiation  from  the  hot  gases  becomes  significant  for 
altitudes  below  25-30  km.^^ 

Figure  4-2  gives  more  details  of  the  flow  over  a  reentry  vehicle.  Figure  4-2(a) 
shows  the  shock  wave  on  the  vehicle,  indicating  the  stagnation  region^^  near  the  blunt, 
rounded  nose  where  the  heating  rate  is  greatest.  To  survive  this  high  heating  rate,  a  reentry 
vehicle  typically  has  a  heat  shield  composed  of  special  material  which  ablates,  i.e., 
decomposes  while  dissipating  the  large  energy  input.  Downstream  from  the  stagnation 
region  the  shock  is  oblique,  i.e.,  weaker,  and  the  boundary  layer  along  the  vehicle 
gradually  thickens;  we  show  schematically  the  transition  from  laminar  to  turbulent  flow  in 
the  boundary  layer,  which  typically  occurs  between  30  and  20  km  altitude.  For  a  typical 
reentry  vehicle  with  an  ablative  heat  shield  the  stagnation  temperamre  is  in  the  range 
2500  K  to  3000  K,  while  the  boundary  layer  temperature  is  in  the  range  1000  K  to  2000  K. 
Figure  4-2(b)  indicates  schematically  how  the  velocity,  temperature,  pressure,  and  density 
vary  as  one  goes  in  the  stagnation  region  through  the  shock  and  to  the  vehicle  wall  "w." 
Note  in  particular  the  structure  in  temperature  through  the  shock.  This  is  determined  by  the 
transfer  of  energy  from  initial  kinetic  energy  into  various  degrees  of  freedom,  which  is 
displayed  in  more  detail  in  Fig.  4-2(c)  which  shows  the  effects  of  kinetics:  translational 
and  rotational  equilibration  occurs  very  rapidly,  giving  a  high  translational  temperature 
immediately  across  the  shock  as  all  the  directed  energy  of  motion  is  transformed  into 
random  translational  and  rotational  energy.  Transfer  into  vibrational  excitation  takes  longer 
and  leads  to  a  lower  temperature  as  the  energy  is  distributed  among  more  degrees  of 
freedom.  Dissociation  and  electronic  excitation  takes  longer  still,  giving  yet  a  further 
reduction  in  kinetic  temperature,  while  ionization  is  slower  still,  yielding  the  final 
equilibrium  temperature.  In  this  sketch,  the  energy  distribution  between  different  degrees 
of  freedom  corresponds  to  a  normal  shock  in  air  at  a  velocity  of  7  km/sec-see  Table  4-1, 
which  shows  that  here  most  of  the  energy  goes  into  dissociation,  a  significant  amount  into 
vibrational  excitation,  and  very  little  into  ionization.^^ 


This  radiation  effect  increases  with  increasing  density  because  the  vehicle  boundary  layer  is  evidently 
optically  thin  down  to  the  lowest  significant  altitudes.  This  whole  issue  will  be  discus^  in  Part  II,  to 
be  prepared  later. 

The  stagnation  region  is  defined  as  that  region  in  the  flow  where  the  ordered  flow  is  brought  to  rest,  so 
that  all  the  ordered  kinetic  energy  of  motion  must  be  transformed  into  thermal  energy.  Thus  it  is  here 
that  the  temperature  and  hence  the  heating  rate  have  their  maximum  values. 

By  contrast,  for  a  nramal  shock  in  Argon  most  of  the  energy  goes  into  ionization. 
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Figure  4-2.  Schematic  Air  Flow  on  a  Reentry  Vehicle. 

(a)  Shock  Wave  and  Boundary  Layer  on  the  Vehicle. 

(b)  Profile  of  temperature,  density,  velocity,  and  pressure  through 
the  Bow  Shock.  (Source:  Murfy,1987-from  OSC). 

(c)  Equilibration  through  Hypersonic  Shock. 
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Table  4-1.  Energy  Balance  Behind  a  Normal  Shock  In  Air  (Source:  Bauer,  1969) 

(A  velocity  of  5  km/sec  corresponds  to  IRBM  entry  while  7  km/sec  corresponds  to  ICBM  entry. 
The  ambient  conditions  are  Po  =  0.01  atm.,  Tq  =  294  K. 
and  energies  are  listed  in  eV/molecule  or  atom.) 


Uo 

5  km/sec 

7  km/sec 

Upstream  of  the  shock 

Directed  kinetic  energy 

Thermal  enthalpy  CpT 

3.76 

0.006 

7.3 

0.006 

Downstream  of  the  shock 

Stagnation  temperature  (K) 

5800  (~  0.5  eV) 

7500  (~  0.65  eV) 

Enthalpy  cf  translation  (3/2)  kT 

1.25 

1.62 

Rotational  enthalpy,  kT 

0.5 

0.65 

Vibrational  enthalpy 

0.37 

0.51 

Dissociation  energy 

1.51 

2.72 

Electronic  energy 

0.011 

0.075 

Ionization  energy 

0.002 

0.006 

Chemical  energy  (N2  +  O2  <=>  2NO) 

-0.02 

-0.01 

Equilibrium  composition  downstream  of  the  shock 

IN2] 

0.579 

0.304 

[NO] 

0.0164 

0.00615 

[02] 

0.000974 

0.000124 

[N] 

0.0840 

0.418 

[0] 

0.319 

0.272 

[e-]  -  [NO-  ] 

0.000137 

0.000630 

6 


4.3  NUCLEAR  FIREBALL 


A  nuclear  weapon  emits  a  significant  portion  of  its  energy  in  the  form  of  X-rays, 
which  are  absorbed  in  the  atmosphere.  The  energy  spectrum  of  these  X-rays  is  normally 
continuous,  having  a  Planck  spectrum  of  temperature^^  1-10  keV.  Figure  4-3  gives  the 
absorption  and  extinction  (=  absorption  +  scattering)  coefficient  of  keV  X-rays  as  a 
function  of  energy  by  "air,"  which  is  four  parts  of  nitrogen  to  one  part  of  oxygen.^^  Note 
that 

(a)  The  extinction  coefficient  falls  off  very  rapidly  with  increasing  energy,  and 
scattering  is  only  important  at  the  highest  energies. 

(b)  A  Planck  blackbody  spectrum  of  temperature  T  has  its  maximum  energy 
density  at  frequency  Vm  such  that^® 

hVm/kT  =  2.825  .  (4.2) 

For  the  extinction  of  X-rays  in  the  atmosphere,  one  can  define  a  dimensionless 
optical  thickness  X 

X  =  (absorption  coefficient)  x  (column  density)  ,  (4.3) 

and  if  X  »1,  the  energy  is  absorbed  in  the  relevant  column.  Table  4-2  gives  the  optical 
thickness  X  for  1  keV  and  10  keV  X-rays  in  a  10-km  horizontal  path  in  air  at  different 
altitudes.  We  see  that  at  altitudes  below  30  km  to  50  km  the  X-ray  energy  is  absorbed  near 
the  burst,  giving  rise  to  a  well-confined  fireball  ( ~  1-10  km),  but  at  altitudes  above  70  km 
to  100  km  the  density  is  so  low  that  energy  deposition  occurs  over  a  very  large  atmospheric 
region  (~  100-1000  km),  giving  rise  to  a  very  different  sort  of  phenomenology. 
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The  notation  is  conventional  but  confusing,  in  that  the  temperature  T  of  a  black-body  spectrum  is 
normally  measured  as  voltage  V  from  the  relation  kT  =  eV,  Note  that  the  differential  energy  spectrum 
of  a  black  body  temperature  1  keV  peaks  at  an  energy  of  2.8  keV-see  Eq.  (4.2) 

Note  that  the  details  of  molecular  structure  are  totally  immaterial  for  incident  X-rays  in  the  keV  energy 
range,  because  chemical  binding  energies-of  5  to  10  eV-are  so  much  smaller  than  the  X-ray  energy. 

The  spectrum  as  a  function  of  wavelength  has  its  maximum  at  X.ni  where  hcAml^T  =  4.965.  TT.is 
gives  the  well-known  Wien's  displacement  law  Xn,  T  =  0.290  cm-K. 
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Source;  Storm  and  Israel  (1970). 


Figure  4-3.  keV  X-Ray  Absorption  and  Extinction  Coefficients  for  Air 
(Source:  Storm  and  Israel,  1970) 


ENERGY  (keV) 


Table  4-2.  Optical  Thickness  for  keV  X-rays  In  a  Horizontal  Path  of  lO-km 
Length  In  the  Atmosphere  at  Different  Altitudes.  (Atmospheric  data  from  USAF, 
1965;  absorption  coefficient  at  1  keV  is  4x10^  cm^/g,  at  10  keV  it  is  3  cm^/g) 
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The  following  description  of  a  low  air  burst  (i.e.,  within  the  atmosphere)  is  taken 
verbatim  from  J.  Zinn  (LANL~see  Barasch,  1979).^^ 

"Within  the  first  microsecond  after  detonation,  the  entire  nuclear  yield  is 
deposited  in  the  bomb  materials,  heating  them  to  extremely  high 
temperatures  in  the  range  10^  K.  Some  of  this  energy  is  immediately 
radiated  away  in  the  form  of  X-rays  and  extreme  ultraviolet  radiation,  but 
this  radiation  is  immediately  reabsorbed  in  the  air  within  the  first  few  meters 
surrounding  the  device.  The  air  is  thereby  heated  to  a  very  high 
temperature,  in  the  range  of  10^  K,  with  a  corresponding  increase  in 
pressure  to  several  thousand  atmospheres.  This  heated  air  and  the  hot  bomb 
vapors  (debris)  in  the  interior  constitute  the  initial  fireball.  The  fireball 
subsequently  expands  through  a  combination  of  radiative  and  hydrodynamic 
processes.  An  intense  shock  wave  (thin  region  of  highly  compressed  air)  is 
formed  at  the  surface,  which  expands  outward  at  a  high  velocity.  Inside  the 
shock,  the  energy  is  rapidly  redistributed  through  emission  and  reabsorption 
of  ultraviolet  rai^ation. 

The  time  variation  of  light  emission  during  the  first  pulse  of  Fig.  4-4  is 
controlled  by  the  radial  growth  of  the  shock.  Although  the  shock  surface  is 
intensely  luminous,  the  shock  is  opaque  to  visible  light,  thus  concealing  the 
air  and  bomb  debris  in  its  interior.  As  the  shock  expands,  it  engulfs  cold 
air,  causing  its  temperature  to  decrease,  with  a  corresponding  decrease  in 
brightness.  TTie  optical  energy  emission  rate  (power)  is  proportional  to  the 
fireball  brightness  times  the  area  of  its  surface.  After  detonation,  the  optical 
power  first  increases,  because  the  increase  in  surface  area  outweighs  the 


^ '  For  a  more  detailed  discussion,  see  Brode,  1%8. 
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decrease  in  brightness.  Later,  it  decreases,  because  the  decrease  in 
brightness  is  dominant. 

During  the  early  expansion  and  cooling  of  the  fireball,  the  ultraviolet 
radiation  in  the  interior  decreases  rapidly  until  it  cannot  effectively 
redistribute  the  energy.  Thereafter,  the  shock  temperature  falls  more  rapidly 
than  the  temperature  in  the  interior.  As  the  expansion  proceeds,  the  shock 
continues  to  cool,  and  as  it  does  so,  it  becomes  less  and  less  opaque  to 
visible  light,  while  also  becoming  less  luminous.  A  point  is  reached, 
corresponding  to  the  minimum  temperature  region  in  Fig.  4-4,  where  the 
shock  becomes  sufficiendy  transparent  to  allow  hght  from  the  hotter  interior 
to  begin  to  escape,  causing  the  optical  output  to  begin  increasing  again. 

Further  growth  of  the  shock  results  in  increasing  transparency  and  a  flmher 
increase  of  the  luminosity  of  the  inner  region. 

As  the  expansion  continues  further,  the  temperature  of  the  inner  fireball 
decreases,  due  to  the  combined  effects  of  hyd^ynamic  expansion  and  loss 
of  energy  by  radiation.  'Phis  cooling  eventually  results  in  the  decrease  in 
luminosity  following  the  second  maximum  in  Fig.  4-4.  After  the  second 
maximum,  the  inner  fireball  gradually  becomes  transparent,  revealing  the 
bomb  debris  inside,  whose  brightness  also  decreases  with  time.  The 
majority  of  the  total  energy  radiated  by  the  fireball  comes  from  the  second 
peak,  which  is  not  very  different  in  its  instantaneous  brighmess  relative  to 
the  first  peak,  but  lasts  about  100  times  longer.  By  the  time  this  "thermal 
pulse"  is  over,  the  available  energy  has  either  been  radiated  away,  trapped  in 
the  debris  and  entrained  air  of  the  fireball,  or  it  is  propagating  in  the  shock 
wave,  which  is  now  well  outside  the  visible  fireball." 

The  fireball  remains  in  the  1000  K  to  10,000  K  temperature  range  for  times  of  3  sec 
to  30  sec  depending  on  bomb  yield,  so  that  its  phenomenology  is  of  considerable  systems 
implication.  For  a  discussion,  see  "The  Effects  of  Nuclear  Weapons"”referenced  as 
ENW,  1964,  Ch.  1  and  2,  in  particular  pp.  87-100,  or  Zeldovich  and  Raizer,  1966,  Vol.II, 
p.  61  If  for  a  description  taken  from  the  first  (1957)  edition  of  ENW.^2 


Zeldovich  &  Raizer,  1966  describe  the  description  of  the  evolution  of  a  fireball  in  the  first  (1957) 
edition  of  ENW  as  better  than  that  in  the  1962/64  edition.  There  is  now  a  third  (1977)  edition  of 
ENW,  which  is  poorer  still  from  the  standpoint  of  describing  what  actually  was  observed  during  a 
nuclear  explosion,  because  by  1977  the  observations  of  atmospheric  nuclear  tests  of  1962  and  earlier 
had  been  massaged  to  provide  a  homogenized  but  not  necessarily  mote  accurate  description  of  what  goes 
on. 
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Variation  of  Apparent  Fireball  Surface  Temperature  with  Time 
In  a  20-kt  Explosion.  (Source:  ENW,  1964,  p.  75) 
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4.4  UPPER  ATMOSPHERE 

Figure  4-5  shows  the  penetration  of  the  solar  radiation  into  the  earth's  atmosphere. 
The  photodissociation  of  O2  is  responsible  for  the  absorption  of  solar  UV  in  the 
150-nm  to  180-nm  (  1500  to  1800  A)  range  which  occurs  in  the  100-km  altitude  range 
and.provides  the  energy  input  which,  in  effect,  drives  the  thermosphere,  while  the 
absorption  of  radiation  in  the  200-300  nm  wavelength  range  by  O3  takes  place  at  altitudes 
near  30  km  to  50  km  and  is  responsible  for  the  stratospheric  photochemistry.  Figures  3-6 
(p.  3-31)  and  3-7  (p.  3-32)  show  the  corresponding  absorption  cross  section  of  the  O2  and 
03  molecules.  Table  4-3  demonstrates  that  all  this  photochemistry  is  driven  by  a  very  small 
fraction  of  the  solar  irradiance;  roughly  half  of  the  total  solar  irradiance  penetrates  the 
atmosphere,  while  only  0.015  of  the  earth’s  atmosphere  is  above  30-km  altitude,  and 
4  X  10“^  above  100  km. 

It  is  rather  easy  to  demonstrate  that  O2  dissociation  in  the  thermosphere  (80  km  to 
300  km,  say)  is  not  in  LTE,  but  rather  in  Radiative  Steady  State  (RSS).  Table  4-4  lists 
temperature,  density,  O2  and  O  densities  and  degree  of  dissociation  defined  as 

<|)diss  =  n(y[no  +  no2]  ('^■'^) 

under  representative  mean  conditions,  and  also  the  LTE  degree  of  dissociation  computed 
from  the  analysis  of  Section  2.2.2.  This  shows  clearly  that  LTE  simply  cannot  explain  the 
degree  of  dissociation  that  is  actually  found, 

A  simple  but  often  useful  way  to  describe  the  non-LTE  upper  atmosphere  without 
going  to  the  full  detail  of  individual  species  concentrations  is  to  introduce  population 
temperatures  for  different  degrees  of  freedom,  as  is  suggested  in  Section  2.5.  Some 
representative  population  temperatures  are  listed  in  Table  4-5.  It  should  be  noted  that- 
because  the  dissociation  and  ionization  occurs  through  the  direct  absorption  of  solar 
radiation— the  dissociation  and  ionization  excitation  temperatures  are  much  larger  than  the 
kinetic  temperatures  of  neutrals  and  electrons,  and  that  different  degrees  of  freedom  (e.g., 
N2  and  O2  dissociation)  have  different  population  temperatures  because  their  excitation 
mechanisms  differ.^374 


^  ^  Missing  data  normally  indicate  that  the  species  are  not  important  at  the  relevant  altitudes. 

For  a  later  version  listing  vibrational  population  temperatures  of  IR  active  molecules-03,  H2O.  CO2, 
CO,  NO— in  the  60-km  to  120-km  altitude  range-see  Degges  and  Smith,  1977,  and  Degges  and 
D'Agati,  1984;  for  a  current  version  of  this  code,  SHARC  (Strategic  High-Altitude  Radiance  Code)  sec 
Sharma  et  al  ,1989. 
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Table  4-3.  Fraction  of  Solar  Irradiance  In  Different  Wavelength  Intervals 

(Source:  USAF,  1965,  p.  16-3) 


Wavelength  Interval 

Fraction  of  Solar  Spectral  Irradiance 

<  121.6  nm 

4.4  X  10~® 

121.6  nm  (H  Ly-a  line) 

3.1  X  10-® 

121.6-180  nm 

2.3x10-5 

180-225  nm 

6.5  X  10-^ 

225-300  nm 

0.013 

300-400  nm 

0.079 

400-700  nm 

0.40 

7000-1000  nm 

0.22 

1-2  tun 

0.23 

>  2um 

0.059 

Table  4-4.  02  Dissociation  In  the  Thermosphere:  Actual  and  LTE. 
(Source  of  input  data:  Bortner  and  Kummler,  1969,  and 
U.S.  Standard  Atmosphere,  1962) 


h  (km) 

n(cm-3) 

T(K) 

n(02)(crTr3) 

n(0)(cmr3) 

<t>diss®°1 

Kdiss(T)* 

80 

3.8  X  10^^ 

181 

8.0  X  10^3 

8.5  X  1010“ 

1.1  X  10-3 

1.6  X  10-117 

4.5  X  10-66 

100 

1.1  X  10^3 

210 

2.0  X  10^2 

5.0x  IOII 

0.20 

6.0  X  10-98 

1.6  X  10-65 

140 

9.5  X  10^° 

714 

9.8  X  10® 

2.2x  IOIO 

0.70 

1.4x  10-11 

8.1  X  10-12 

200 

6.9  X  10^ 

1236 

3.6  X  108 

3.2  X  10® 

0.90 

2.5  X  10^ 

3.4  X  10-3 

300 

6.7  X  108 

1432 

1,1  X  10^ 

5.1  X  108 

0.98 

1.8  X  107 

0.13 

*  K,jiM(T)  is  the  O2  dissociation  equilibrium  constant  of  Eq.  (2.2). 

**  At  80  km  and  below,  the  O-atom  number  density  at  night  is  less  than  in  the  day  time,  by  a  large  factor. 
Daytime  values  are  listed  here. 
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Table  4>5.  Population  Temperatures  for  Different  Degrees  of  Freedom  In  the 
Upper  Atmosphere  (Source:  Bauer,  Kummler  and  Bortner,  1971) 


Altitude  (km) 

100 

200 

300 

Kinetic  temperature  of  neutrais  (K) 

185 

930 

1090 

Kinetic  temperature  of  electrons  (K) 

370 

1570 

1900 

O2  dissociation  temperature  (K) 

1900 

1870 

1830 

N2  dissociation  temperature  (K) 

1800 

2800 

(2400) 

NO  ionization  temperature  (K) 

2560 

2850 

- 

O2  ionization  temperature  (K) 

2630 

3100 

- 

0(^D)  excitation  temperature  (K) 

1280 

1870 

- 

O2  (a ‘'4)  excitation  temperature  (K) 

1000 

- 

- 

The  physics  of  the  upper  atmosphere  will  not  be  discussed  in  detail  here.  A 
discussion  emphasizing  the  overall  physics  and  simple  but  usable  numbers  is  given  in 
Bauer,  1972.  There  exist  much  more  modem  (and  authoritative)  treatments,  our  current 
understanding  of  the  stratosphere  is  much  better  than  it  was  in  1972  and  also  my  discussion 
of  meteorology  as  given  there  is  bad.  However,  while  we  have  much  more  data  on  the 
thermosphere— see  e.g.  U.S.  National  Report  to  lUGG,  1987— yet  our  understanding  of  the 
basic  physics  has  changed  very  little,  and  nxxlem  books  have  a  tendency  to  overwhelm  the 
non-expert  reader  with  details. 
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This  appendix  presents  a  brief  treatment  of  the  quantum  mechanical/spectroscopic 
structure  of  various  air  species. 

A.l  ATOMS 

The  electronic  structure  of  atoms  is  built  up  by  filling  electron  shells  subject  to  the 
Pauli  exclusion  principle,  which  says  that  for  given  principal  quantum  number  n  and  orbital 
angular  momentum  quantum  number  1  (where  1  <  n)  one  can  have  2(21  +  1)  electrons.^  The 
notation  is  odd  but  traditional^: 

total  orbital  angular  momentum  (1  or  L)  0  1  2  3  4  5 

notation  s(or  S)  p(or  P)  d(or  D)  f(or  F)  g(or  G)  ... 

where  (nl)  characterizes  an  individual  electron  wave  function,  while  L  is  the  orbital  angular 
momentum  of  the  atom.  In  other  words,  for  given  principal  quantum  number  n  we  can 
have  2  s-electrons  (1  =  0),  6  p-electrons  (1  =  1),  10  d-electrons  (1  =  2),  etc. 

Thus,  the  atoms  in  the  periodic  table  are  built  up  like  this: 


Z  =  1 

H: 

(Is)  2s  (where  2  refers  to  a  doublet,  i.e.,  one  electron  of  spin 
s  =  i/2,  2s  +  1  =  2) 

II 

N 

He: 

(ls)2  Is  ("singlet  S,"  ^S,  is  a  closed  shell,  two  electrons  of 
opposite  spins) 

^  Each  electron  has  four  quantum  numbers,  which  may  be  taken  as  the  principal  quantum  number  n,  the 
orbital  angular  momentum  quantum  number  I,  the  (integral)  component  Iz  of  1  ^ong  some  designated 
direction  z,  and  the  spin  quantum  number  Sz  corresponding  to  angular  momentum  +/-  h/2jt  along  the 
designated  z-direction.  The  exclusion  principle  says  that  no  two  electrons  in  an  atom  can  have  all  four 
of  (n,Uz.Sz)  ihe  same. 

2  S  stands  for  "sharp,”  P  for  "principal,"  D  for  "diffuse,"  and  F  for  "fundamental,"  empirical  terms 
applied  to  the  characters  of  spectral  lines  in  the  late  1890s  as  spectrographs  were  first  used 
systematically  to  study  atomic  and  molecular  spectra. 
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Z  =  3 


U 


I 

Z  =  4  Be 


Z  =  5  B 

Z  =  6  C 

Z  =  7  N 


Z  =  8  O 


Z  =  9  F 


Z  =  10  Ne 
Z  =  ll  Na 


(1s)2(2s)  [one  (2s)  electron  outside  a  closed  shell,  i.e.,  doublet 
or  unpaired  electron] 

(1s)2(2s)2  [two  (2s)  electrons  of  opposite  spin,  (2p)  comes  next, 
room  for  6  of  them] 

(ls)2(2s)2(2p)  2p 

(ls)2(2s)2(2p)2  3p  ("triplet  P"— the  two  electron  spins  are  parallel, 
2s  +  1  =  3  if  s  =  1) 

(ls)2(2s)2(2p)3  ("quartet  S,"  s  =  3/2.  From  (2p)3  one  can  also 
make  2d  and  2p  configurations,  which  are,  respectively,  2.38  eV 
and  3.58  eV  above  the  ground  state.  Figure  A-1,  which  shows 
the  Potential  Energy  curves  for  N2,  and  Fig.  A-2,  which  shows 
those  for  NO,  indicate  molecular  states  involving  N(2D)  and  N(2p), 
as  weU  as  t  those  that  dissociate  into  N('^S)) 

(ls)2(2s)2(2p)^  3p  (from  (2p)‘*  one  can  also  make  and 
configurations,  which  are,  respectively,  1.97  and  4.19  eV  above  the 
3p  ground  state,  where  the  triplet  splitting  is  0.02  and  0.04  eV. 
Figures  A-2  for  NO  and  A-3  for  O2  indicate  molecular  states 
involving  0(^D)  and  0(1  S).  They  also  show  molecular  negative 
ions  dissociating  into  the  atomic  ion  O"  (2po)  which  has  9  electrons 
outside  the  nucleus  and  is  thus  isoelectronic  with  Ruorine. 

(ls)2(2s)2(2p)5  2p  (one  electron  missing  before  a  closed  shell 
(halogen) 

(ls)2(2s)2(2p)^  1$  (closed  shell-noble  gas) 

(ls)2(2s)2(2p)^(3s)  2s  (next  alkali,  one  electron  outside  a  closed 
shell). 


etc. 

Atoms  have  electronically  excited  states,  such  as  the  Rydberg  states  of  the  H-atom, 
and  of  other  atoms  having  just  one  electron  outside  a  closed  shell,  but  these  states  lie  10  eV 
or  more  above  the  ground  state,  so  that  transitions  involving  them  lie  in  the  far  UV  (10  eV 
corresponds  to  1240  A  or  124  nm  wavelength).  For  viewing  through  the  atmosphere  we 
are  mainly  interested  in  near- visible  radiation,  i.e.,  X  ~  300  nm  to  700  nm  or  energies  of 
4.1  eV  to  1.8  eV.  Important  transitions  in  atomic  oxygen  correspond  to  the  dipole- 
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forbidden  630  ran  "red"  (^D  =>  3p)  and  558  nm  "green"  (^S  =>  ^D)  auroral  lines,  and  to 
the  atomic  nitrogen  2p  =>  2d  (1040  nm)  and  2d  =>  (520  nm)  lines.  Figures  A-1  to  A-3 

show  that  the  air  molecules  N2,  NO,  and  O2  have  many  low-lying  states  with  important 
transitions  between  them,  so  let  us  now  talk  about  molecular  transitions. 


t‘ol«nl)al  •n«/gy  laltcifon  vollt) 


Poltnlial  anargy  (•Icclron  volli) 


0.4  0.8  1.2  1.6  2.0  2.4  2.8  3.2  3.6 

Initrnuclaar  distonct  (A) 


•  Figure  A-2.  Potential  Energy  Curves  for  NO  (Source:  Gilmore,  1964) 
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Polanlioi  ana/yy  (tiiciioii  «oiiSj 


Figure  A-3.  Potential  Energy  Curves  for  O2  (Source:  Gilmore,  1964) 


A. 2  DIATOMIC  MOLECULES 


Here  a  key  parameter  is  the  angular  momentum  along  the  molecular  axis:^ 
angular  momentum:  0  12  3 

notation:  Z  O  A  O  . 

Another  key  parameter  is  the  spin  s,  or  equivalent  multiplicity  (2s+l);  a  molecular  state  is 
labeled  A.  The  (electronic)  ground  state  of  a  molecule  is  always  called  "X,"  with 
excited  states  labeled  A,B.—  (in  increasing  energy^).  The  important  air  molecules  N2, 
NO,  and  O2  have  very  complex  sets  of  electronically  excited  states;^  which  are  shown  in 
Figs.  A-1  to  A-3  respectively. 

H2:  H(ls  ^S)  H(ls  ^S)  gives  a  ground  state  X'Z. 

N2:  Ground  state  is  X^Zg'*'. 

Figure  A-1  (from  Gilmore,  1964)  suggests  that  A^Zy^  is  the  only  other  stable  state  that 
dissociates  into  N(4S)  +  NC^S).  However,  current  calculations  and  observations  (Partridge 
et  al.,  1988)  suggest  that  the  5Zg+  state  is  bound  more  strongly  than  suggested  in  Fig.  A-1, 
and  that  in  fact  this  A"  state  is  the  primary  precursor  state  involved  in  the  yellow  Lewis- 
Rayleigh  afterglow  of  N2. 

NO:  Ground  state  X^O,  See  Fig.  2-6  [a'^FIj  also  dissociates  into  N('^S)  +  0(3p)] 
O2:  Ground  state  X^Zg”.  See  Fig.  2-7  . 

The  other  indexes  correspond  to  different  symmetries  of  the  atom  or  molecule,  which  are 
important  for  radiative  selection  rules.  In  the  case  of  an  atom,  if  Zjli  is  odd  or  even  the 
state  is  described  as  odd  or  even,  respectively,  while  for  a  homonuclear  diatomic  molecule 
where  the  electronic  wave  function  either  remains  the  same  or  changes  sign  when  it  is 
reflected  at  the  center  of  the  interatomic  line  we  describe  the  state  as  even  (g)  or  odd  (u).^ 


^  The  convention  is  the  following:  s,p,d,...  applies  to  individual  atomic  electrons;  S,P,D  applies  to 
(multi-electron)  configurations  in  an  atom;  a,Ji,5  applies  to  single-electron  states  in  a  molecule,  and 
I,n,4  applies  to  (multi-electron)  configurations  in  a  molecule. 

^  Excited  states  of  the  same  multiplicity  (spin)  as  the  ground  state  are  labeled  A3,  etc.,  while  those  of  a 
different  multiplicity  are  labeled  a,b,  etc.  N2  is  an  exception;  the  convention  is  reversed  due  to  a 
historical  error. 

^  Principally  because  the  ground  states  of  the  N  and  O  atoms  are  not  but  rather  and  ^P, 
respectively. 

^  This  property  is  called  parity. 
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Further,  in  the  case  of  a  symmetric  top  (homonuclear  diatomic  molecule)  there  are  two 
wave  functions  of  different  symmetry,  labeled  and  respectively.’^ 

As  an  example,  reference  to  Fig.  A-3  shows  that  two  0(^P)  atoms  can  come 
together  along  different  potential  energy  curves  to  give  a  variety  of  different  molecular 
states,  both  stable  ones  of  different  spin  such  as  X^Lg",  C^Aq  and  a^Ag,  blLg+,  and  also 
unstable  (repulsive)  states  such  as  ^>3  Ilu  g,  etc.  The  important  "Schumann-Runge" 
transition  (discussed  in  Chapter  III)  is  a  transition  =>  X^Eg”  which  demonstrates 

several  selection  rules: 

(i)  spin  remains  constant  (triplet  to  triplet) 

(ii)  parity  changes  (u  =>  g) 

(iii)  inversion  symmetry  remains  constant  (-  =>  -). 

A. 3  TRIATOMIC  MOLECULES 

Things  get  more  complicated.  A  triatomic  molecule  may  be  linear,  like  the  CO2 
ground  state  (^E),  or  bent,  e.g.,  NO2,  O3,  H2O.  The  ground  and  excited  states  are 
designated  by  ~X',  "A',  ~B',  etc.  For  linear  molecules  the  same  state  classification 
symbols  (L,  IT,  etc.)  are  used  as  for  diatomic  molecules,  while  for  nonlinear  molecules  the 
electronic  symmetry  is  designated  by  the  symbols  Ai,  A2,  Bi,  and  B2  for  symmetric 
molecules  like  H2O  and  NO2,  and  the  symbols  A'  and  A”  for  nonsymmetric  molecules  like 
HNO.  A  non-linear  triatomic  molecule  has  three  fundamental  vibration  frequencies,®  and 
different  electronic  states  may  be  straight  or  bent  differently  from  the  grour*  state.  For  a 
discussion  of  the  structure  of  these  molecules,  see  Herzberg,  1966;  for  the  applications 
discussed  here  it  is  sufficient  to  know  that  there  are  different  types  of  states  (electronic, 
vibrational,  and  rotational,  with  relatively  weak  coupling  between  them). 


^  Note  that  the  statistical  weight  of  a  state  is  15,  while  that  of  a  state  is  6.  In  each  case  the 
superscript  3  corresponds  to  spin  1,  i.e.,  multiplicity  3,  but  D  corresponds  to  L  =  2,  i.e.,  statistical 
weight  (2L+1)  =  5.  By  contrast.  A  (A.>0)  corresponds  to  a  symmetric  top  with  doubly  degenerate 
rotation^  levels-see,  e.g.,  Herzberg,  1950,  p.  128. 

®  See  footnote  17  on  p.  3-1  of  the  main  text. 
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For  the  dissociation  of  a  diatomic  molecule,  one  may  use  the  standard  results— see 
any  text  on  statistical  mechanics,  e.g.,  Fowler  and  Guggenheim,  1949,  pp.  87, 96,  165: 

f(J)  =(27iM(J)kT/h2)3/2vgjrj  (J  =  A,B)  (B.la) 

f(AB)  =  (27tM(AB)kT/h2)3/2  v  r(D  q(T)  gAB  ta  rs  /  Cab  (B.  lb) 

r(T)  =  rotational  partition  function^  =  T/0r  ;  ©r  =  hcBo/k  (B.  Ic) 

q(T)  =  vibrational  partition  function2  =  l/[l-exp(-©v/T)];  ©y  =  hccOo/k,  (B.  Id) 

and 

gj  =  electronic  partition  function  (statistical  weight ) 
rj  =  nuclear  spin  weight^ 

<^AB  =  symmetry  number  =  2  if  A=B,  otherwise,  1. 

Bo,a)o,To,  etc.,  come  from  Huber  and  Herzberg,  1979,  or  equivalent  and  are 
measured  in  cm“^. 

Let  us  do  calculations  for  the  O2  molecule  for  which: 

D  =  5. 1 16  eV,  i.e.,  ©d  =  59300  K 

(Do  =  1568  cm-l,  i.e.,  ©v  =  (hcA)cOe  =  1.439  cOq  =  2256  K 

Bo  =  1 .438  cm-1,  ©r  =  (hc/k)Bo  =  2.07  K 


^  This  is  an  approximation  good  when  the  rotationaJ  level  spacing  is  small  compared  to  kT,  which  holds 
for  molecules  such  as  N2, 02.  NO,  etc.  (not  H2!)  at  temperatures  above  ~  200  K. 

2  An  approximation  good  for  harmonic  oscillators  with  vibrational  spacing  large  compared  with  kT. 

3  For  the  dissociation  of  a  diatomic  molecule,  the  term  rAB/fArs  cQual  unhy  because  the  nuclear 
spins  do  not  change-see,  e.g.,  Fowler  and  Guggenheim,  1949,  p.  166. 
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g02(T)  =3  +  2  exp(-©E/T) ;  €>e  =  (hc/k)Tei  =  1 1400  K  (contribution  from 
lA  state) 

go  =  3  . 

Calculations  for  various  values  of  number  density  ratios  n/no  .where  no  =  2.687  x 
cm~3,  respectively)  are  shown  in  Fig.  2-1,  where  we  list  a  fractional  dissociation  defined 
as 

<l><liss  =  JVyfno  + 1^02]  •  (B.2) 
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APPENDIX  C 

EINSTEIN  A-  AND  B-COEFFICIENTS,  f-NUMBERS  AND 

BAND  STRENGTHS 


This  appendix  presents  the  relation  between  the  various  measures  of  band  strength, 
following  Einstein's  derivation  of  the  A  (spontaneous)  and  B  (stimulated)  emission/ 
absorption  coefficients  as  presented  by  Fenner,  1959-see  also  Ludwig  et  al,  1973.  We 
use  cgs  units  here  because  this  is  conventional. 

Consider  a  transition  from  an  upper  state  u  to  a  lower  state  1,  with  Nj  molecules  in 
state  j  (=  u,l)  per  unit  volume.  Ayi  and  Bui  are  emission  coefficients  (spontaneous  and 
stimulated)  and  Biu  is  the  stimulated  absorption  coefficient.  The  spectral  density  of  a  black- 
body  radiation  field  is:^ 


p(v)  =  (87thv3/c3)  (ehvAT ..  i)-l  ^ 

(C.l) 

whose  units  are  erg  sec/cm^,  while  Bjj  is  measured  in  cm^/erg  sec^ 
measured  in  sec“^ 

=  cm/gm  and  Ajj  is 

In  Radiative  Steady  State,  defined  as 

Absorption  =  Emission  , 

(C.2) 

we  have 

Ni  Bju  p(vju)  =  Ny  [  Aui  +  By]  p(viu)] 

(C.3) 

and 

Nl/Ny  =  [Ayl  +  Byl  p]  /Bly  p  =  (gl/gyle^V/kT  ^ 

(C.4) 

which  gives  the  following  relations  between  the  A-  and  B-coefficients: 

Ayl  =  (87Chv3/c3)  Byl  , 

(C.5) 

gl  Blu  =  gy  Byl  . 

(C.6) 

'  In  the  notation  of  Section  3.2  of  the  main  text,  one  should  write  B(v,T)  instead  of  p(v),  but  the  use  of 
p(v)  is  conventional  here,  and  there  would  be  confusion  between  blackbody  radiance  and  the  Einstein 
B-coefncient 
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As  in  Eq.(3.34a)  in  the  main  text,  we  define  a  band  strength 
S'iu=Jkv[l-€"''''^’^dv  =  (c2/87cviu2)NiAui(gu/gi)[l-exp(-hviu/kT)]  .  (C.7) 

The  factor  [l-exp(-hviu/kT)]  is  a  correction  for  stimulated  emission,  which  arises  because 
of  the  B-terms  (stimulated  emission  and  absorption).  This  term  is  often  close  to  one  at 
room  temperature  and  can  frequently  be  neglected,  but  it  can  be  important  here.  The 
dimensions  of  band  strength  are 

S'iq  =  (number  density)  x  (absorption  cross  section)  (C.8) 

and  it  is  normally  measured  either  in  10-20  cm- 1 /(molecules  x  cm-2)  or  in  (cm-2  atm-1), 
where... 

S(cm-2  atm*^)  =  2.687  x  (T</r)  S  (cm-Vmol.cm-2)  (C.9) 

with  To  -  273.2  K-  see,  e.g.  Raud  &  Camy-Peyret,  1975,  p.  284.  One  may  write 

Siu(cm-2  atm-1)  =  3.21  x  1028  Aui/Vui^  ,  (C.  10) 

where  Aui  is  measured  in  sec”!  and  Vui  in  Hz. 

One  may  also  define  a  dimensionless  absorption  oscillator  strength  fiu  =  fabs  by 

S'lu  =  (7te2/mc)  Ni  fiu  [  1-exp-hviu/kT]  ,  (C.  1 1 ) 

or  at  Standard  Temperature  and  Pressure  (STP),  if  practically  all  the  absorbing  molecules 
are  in  the  ground  state  1,  provided  that  exp-hviu/kT  «  1 , 

Siu  (cm-2  atm-l)  =  2.3795  x  lO?  fiu  ,  (C.  12) 

and  an  emission  oscillator  strength  fui  is  related  to  the  absorption  oscillator  strength  fiu  by 
the  relation 

fuI  =  (gl/gu)flu  •  (C.  1 3) 

From  (C.IO)  and  (C.12)  we  get 

Aui  =  7.413  X  10-22  Viu2  flu  .  (C.14) 

Finally,  the  lifetime,  Xu,  of  an  excited  state  u  is  related  to  the  A-coefficients  of  all 
lower-lying  states  by  the  relation 

1/Xu  =  2i  Aui  +  Zi  Bui  p(Viu)  +  21u' Buu'P(Vuu')  (C.  15) 
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appendix  d 

EXCERPTS  FROM  PHYSICAL  CONSTANTS,  USAF  1985 

AND 

U.S.  STANDARD  ATMOSPHERE,  1962 


APPENDIX  D 

EXCERPTS  FROM  PHYSICAL  CONSTANTS,  USAF  1985 

AND 

U.S.  STANDARD  ATMOSPHERE,  1962 


Physical  Constants  (from  USAF,  1985,  p.  A-T) . D-4 

U.S.  Standard  Atmosphere,  1962  ("USSA-62"**) . D-5 


*  Agree  to  four  figures  or  belter  with  Cohen  and  Taylor  (1987),  which  is  the  latest  compilation  of 
physical  constants. 

**  Note  that  USSA-62  corresponds  to  near-solar  maximum  conditions  (the  largest  ever  was  in  1957,  with 
Zurich  Sunspot  number  190);  USSA-76  corresponds  to  near-solar  minimum  conditions,  with  a  density 
smaller  by  perhaps  a  factor  3  at  500  km  altitude.  CIRA-72  has  intermediate  values.  For  a  current 
listing  of  reference  and  standard  atmospheres,  see  Whitten  &  Vaughan,  1989. 
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y\ 

UNITS,  CONSTANTS,  AND  CONVERSION  FACTORS  | 


Qijjptity  Symbol  (expression)  Value  in  Si  (cgs)  units  Error  (ppm) 


1  speed  0^  iigrtt  tn  vacjum 

c 

2.997  924  58 X  10*  m-s  '(I0'°cm  sec  ’) 

0.004 

2  Elementary  cnarge 

e 

1  602  189  2x10  ’»C(10  “emu) 

2.9 

(4.803  242x10  '“esu) 

2.9 

3  PlancK's  constant 

n 

6.626  176x10  **  J  s  (10  erg  sec) 

5,4 

fl  =  /l/2ir 

1,054  588  7x10  “J  s(iO  erg  sec) 

5.4 

4  Electron  rest  mass 

m. 

0.910  953  4x10  “kg  (10  "  gm) 

J.t 

S  Avogadro  constant 

N. 

6.022045  X  10“  mol  '  (10“  mol  ') 

5.1 

recent  value 

6,022  097  8  X  1 0“  mol  -  ’  (1  O’*  mol  ' ) 

1.0 

6  Molar  gas  constant 

R 

8.314  41  X  10“  J  mol  '  K  ‘  ’  (10' erg  mol  '■’K  ') 

31 

7  Boltzmann  constant 

k  =  R/N^ 

1.380  662x  10  ”  J-K-'  (10  '*  erg  'K  ') 

32 

a  Gravitational  constant 

G 

6.672 Ox  10  "  N  m’  kg  ’  (10  *  dyn  cm’  gm  ’) 

615 

9  Molar  volume  (273.15 'K, 

31 

Po  =  t  atm) 

k',.  =  ^  ViFo 

22.413  83x  10  *  m*  mol '  ’  (10*  cm’  mol ' ') 

10  Faraday  constant 

F  =  N,e 

9,648  456x10*  C  mol  '  (10*  emu-mol ' ') 

2.8 

1 1  Rydberg  constant 

“(Airfol  ’(m.eV4»ypc) 

1,097  373  177x10' m  ’  (10*  cm  ’) 

0.07 

recent  value 

1.097  373  147  6x  10' m  ’  (10*  cm  ') 

0  0003 

12  Fine  structure  constant 

a  '  =  (4Mo](6c/e^) 

137.036  04 

Oil 

recent  value 

137.035  963 

13.  Classical  electron  radius 

f.  =  [‘‘irroi 

2.817  938  Ox  10  '*  m  (10  ’*  cm) 

2.5 

14  Specific  electron  charge 

e/m. 

1,758  804  7 X  10”  C  kg  '  (10' emu  gm  ’) 

2.8 

15  Electron  Compton  wavelength 

A,  =  A/m,c  =  a  V, 

3,861  590  5x10  ’*  m  (10  ”  cm) 

1.6 

16  Bohr  radius 

ev 

li 

o 

0.529  177 06 X 10 -’®m  (lO  'cm) 

0.82 

1 7  Magnetic  flux  quantum 

=  Id  '(/ic/2e) 

2.067  850  6x10  '*  T  m*  (10  '  Gsem’) 

2.6 

h/e 

4  135  701 X 10' ’*  J  s-C'’  (10"' erg-sec  emu  ') 

2.6 

18.  Quantum  of  circulation 

h/2m. 

3.636  945  5  X  1 0  ■  *  J  Hz  ■ '  kg  '  (1 0“  erg-sec  gm 

') 

1,6 

19  Atomic  mass  unit 

1  u  «  gm-mol  V/V, 

1.660  565  5X10  ”'  kg  (10  “  gm) 

5.1 

20  Proton  rest  mass 

m. 

1.672  648  5x10- "  kg  (10  “  gm) 

5.1 

1  007  276  4  70  u  (amu) 

0.011 

1836  151  52 

0,38 

21.  Neutron  rest  mass 

'T’r, 

1.674  954  3x10  ■*' kg  (10  “  gm) 

5.1 

1.008  665  012  u  (amu) 

0.037 

22.  Electron  g  factor 

iff. 

1.001  159  656  7 

0.0035 

recent  value 

1  001  159  652  200 

0.0004 

23  Bohr  magneton 

fig  =  [c]{e*/2w,c) 

9.274  078x10  “  J-T  '  (10  ”  erg  Gs  ’) 

3.9 

24  Nuclear  magneton 

fi^  =  [cl(ei»/2m,c) 

5.050  824x10-*'  J-T  '  (10  **  erg-Gs  ’) 

3.9 

25  Electron  magnetic  moment 

9.284  832x10  “J-T  '  (10-*’  erg-Gs  ') 

3.9 

26  Proton  magnetic  moment 

fig 

1.410  617  1  X  10-**  J-T  ■ '  (10  **  erg-Gs  ') 

3.9 

p./p. 

658.210  688  0 

0.010 

27.  Proton  gyromagnetic  ratio 

rp 

2.675  198  7  X  10*  rad  s  '-T  '  (10*  rad  sec  '  Gs 

') 

2  8 

26  Stefan-Boltzmann  constant 

a  =  (ir^/60)/C‘/flV 

5.670  32x10-*  W  m  -*  K-*  (10-*  erg  sec  '  em 

*--K  *) 

125 

29  First  radiation  constant 

C,  w  Zrrhc‘ 

3.741  832x10  '*  W  m*  (10 -*  erg-sec  ’-cm*) 

5.4 

30  Second  radiation  constant 

c,  =  hc/k 

1.438  786  X  10  *  m-K  (10°  cm-’K) 

31 

Energy  equivalants 

Quantity 

Symbol  (expression) 

Value 

Emor  (ppm) 

Atomic  mass  unit 

u 

931.501  6  MeV 

2.8 

Proton  mass 

m. 

938  279  6  MeV 

2.8 

Neutron  mass 

m. 

939,573  1  MeV 

2.8 

Electron  mass 

Pf. 

0.511  003  4  MeV 

2.8 

Electron  volt 

1  eV/* 

11  604  50  K 

31 

1  eV/nc 

8  065  479  cm '  ’ 

2.6 

1  ev/n 

2.417  969  6x10’*  Hz 

2.6 

1  ev 

1.602  189  2x  10  ”  ergs 

2.9 

Planck's  constant 

fl 

6.582  173x10  **  MeV  sec 

2.6 

Ac 

1.973  285  8x10  ”  MeV  em 

2.6 

(Ac)’ 

0.389  385  7  GeV*.mb 

52 

Rydberg  constant 

fl.  he 

13.605  804  eV 

2.6 

Voltage-wavelength  product 

vx. 

12  398.520  eV  A 

2.6 

Gas  constant 

R 

1.987  19  cal-mol  ’  'K  ’ 

31 

D-4 


litlluwliiK  ciiliy  kliuultl  Im:  imiUi|tlictl.  (rmillMuctll 


(CtfNllMMrii) 


3  J  “  y  rr 


d?-  * 


Cw^S'S  l^»a-'4-^  ‘Si****<*^f^  #“«»«CS-» 

s  —  —  5af«»«i  ^4?.  rt23irt?<  -H  p*—-?!** 

3*«Jt  -j-j^w^s/?  i/5  —  ssssxr^ 

3<4r4^4^«^4  r474M74  ^^4^94^  ^  C*  94  7^  94  S4  74  94  C4  ^  94  94  94  ^ 


Si**  •^  —  **a*»  a«  —  —  — ••«-•-•»>«•  'iSsSS'-* 

••  —  *•  P»  —  •M*''^  5  —  •4f«-r»  i**u*s:^  -rft  —  —  wiS  S.-??a^^  rty2—  -» 

»  ••?«•^4C^«  •••?—  p*»wt^  —  p-  — >“PS'^C  p»w^  —  2  —  ■M-»c  S'  —  —  a» 

^  -ip's-.  —  •  —  —  -^5*4  94  ^  >A  '■C  e4i'?4r4Jn4P  'f'irt-jiK*  — 1^  —  ^i*>4«'4ei9 


o  —  ^ 

-r  I  I 


—  M  ey 

I  I  I 


^  2  i  1^  *4 


—  sSnS^^i—  —  s»*isap4  a<.$i«  —  9 

p-snmS  p*s>*-3yS  n  ^  94  ^  «4)r»3<^2  m  —  syr-4  t-5N?4^-* 

^Wtwl^  —  ^  ^  si  9^  4i  ^  ^  ^  9i  94  94  94  mt  mi  mi  ^  S^  X  p^i 


?  S'  C*  S'  s 


p*'daas  ai<*p*<^a  — 44e4«a  M«/)p->aa  s^si'I'^o  «*ivS9e^4 

S>pS  —  —  •  w*fa  —  *^-S'  p%p*^»<p*  p«*’(^i^'i<,jf  —  —1^^  —  —  9i/9P*d 

2ss—  s  —  —  wn  nSp»*ss  —  ys— I— 5s  s5  — 

SSiS  ss2s  —  —  ««»••«  «l  'MP4^4SIM 


=  =  »5i 

■c  —  a  vt  *! 

M  S4  —  —  — 


22^2  ?  S  —  2ss 

—  a2a  a  cp»p»x—  p»  —  .5  —  S'  r»r«»p»}'?i 

—  —•XC-  >«  Pt?-SI  ^1  —  —  —  sv  X<C,/i<»^^  P1«nMr4»  ^  mi  mi  si  Xri'Swiwi 


y?  —  2  -•  S  p* 


S'  —  ‘S'rtS  'Sr»Mx 

y:  —  ^  —  ^  ^  —  5\  — 

—  -caxs  XP-P.P* 

—  S'X'S'A  y9^w«») 


I  11  I  III  I 

X  — \rtX  f»— yi  —  «j^  'CSP^Xy'  P'»S*— Sx  X'SX**!  wie*  —  —  p* 

*n  f".  —  •!  —  y*a  p«*r-p»^s;  p*  —  —  *?  X  —  ^5^4  — -csxp^ 

wi  tAsic'  s'^xp*^  p«cs2  —  y?x  —  2  —  p*piS3<rt 

Cy  mi  ^  si  XXWlVpi  94  94  9*  mi  mi  ^  ^  pi  ^  94C4e'i  —  — 


^  X  N*  ^  -C 
••  —  S4  p-  -r 

s  X  s  ••  y* 
x^p-y*  — 
«i  PS  M  si  ?« 


53??^ 


ss;^Ss 

X  s  >4  2! 

«p  X  p^  y;  X 
(*•  P»  X  X  9 


^  P4  pS  tfi 

S  si  ?  s 


••  yiwt^  e»ysN5  ys9>sxsi  — p*»xp»wspi»  x*r  —  yip» 

PI  ysr«—  >s«pPSy3X  e49s4C«A  'Xys— p«s4  S'SX'^pi  yss4nss 

2?  *  '>5  *•  «  ^  vs  X  '4  ^  ^  ^  9  iff  9  y?  9’  ••  p?  p?  ••  ••  y»  4  x  X'  si  >4 

^  ^  ^  9  C  —  **'*pvs  ’3p»p«xCv  ^x25  —  -•  —  ——  -• 

S'  XXO  ^  ^  »m  tm  •m  ^  ^  ^  ^  e4P4S4«VS4 


evfvix  — 

94  94  .A  s  S 


t\  sil5^ 


s’  S  S 

—  ^  f#  fV 


OSXXX  —  9^X—  ©yj  —  ^P»  0>s«xx  S4»fifvO?'4  P>XXXS# 

•PXP^yr4  xyivsp*—  •?•  —  —  S40  P-^PSXP?  —  CsiP**rt  p«33p  —  P4 

*<3^  —  sp»  pis'^xsl  x'xsixsr  xpsxrip-i  p>^p«ip«^p^9  «r*  —  ?p4s 

4S‘SJp»^  p»|sx  —  PS  5— y»ys4  3p*p-XX  XXX9S  xsxsJ' 


I9SS3 

3  s  S 

pw  p»  a 


§333  333i3  §3333  33|i3  §3331  §3111  §333^ 
is§|  sama  SSSlI  islll  |  =  |5i  isiil  |ai3Si 

S*S4e4C«C4  CVCVSS^^  XPSPSPSP)  «PSP1PS^ 


D-7 


